AD-A248  081 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


THESIS 


^DT!C 

,  30. 1992  f 


evaluation  and  improvement 

OF 

THE  ASW  SYSTEM  EVALUATION  TOOL 
by 

Peng-tso  Chang 
March  1992 


Thesis  Advisor: 
Co-Advisor 


Yuh-jeng  Lee 
James  N.  Eagle 


Approved  for  public  relensc;  di^fnhntinn  i.s  unlimited 


92-07837 


UNCLASSIFIED 


SECURmr  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 

1*.  REPORTSECURlTYCLASSinCADON 

UNCLASSIFIED 

lb.  RESTRlCnVEMARKINGS 

2m.  security  CLASSIFICADON  AUTHORITY 

3.  DISTRIBUTION/ AVAILABIUTY  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited. 

2b.  DCLASSIHCADON7DOWNGRADING  SCHEDULE 

.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


3.  MONITORING  ORGANIZATION  REPORT  NUMBER® 


6*.  NAME  OF  PERFORMING  ORGANEAmN  6b.  OFFICE  S^T-IBOL  7«.  NAME  OF  MONITDRING  ORGANIZATION 

Computer  Science  Dept .  W ^pticabU)  Naval  Postgraduate  School 

Naval  Postgraduate  School  CS 


6c  ADDRESS  (ciry,  itaie,  end  2JP  code) 
Monterey,  CA  93943-5000 


8a.  NAME  OF  FUNDING/SPONSORING  6b.  OHICES^'MUOL  9.  PROCUREMENT  INSTRUMENT  IDENTinCATION  NUMBER 

ORGANIZATION'  Of  Applicable) 


8c  (city,  state,  andZlP  code) 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


PROEICr 

NO. 


WORKUNTT 
ACCESSION  NO. 


TITLE  (Include  Security  Classification) 

EVALUATION  AND  IMPROVEMENT  OF  THE  ASW  SYSTEM  EVALUATION  TOOL 


13b.  TIME  COVERED 


IS,  PACE  COUNT 


12.  PERSONAL  AtTHORCS) 

Penq-tso  Ghana 


13*.  TYPE  OF  REPORT 

Master's  Thesis 


16  SUPPLE.MENTARY  NOTATION 

The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the 


14.  DATE  OF  REPORT 
March  1992 


17. 

COSAD  CODES 

FIllD 

GROUP 

SUBGROUP 

18.  SUBIECT  TERMS  (continue  on  reverse  if  necessary  and  identify  by  block  number) 

Object-oriented  Programming,  Ob ject-oirented  Simulation, 
ASW  System  Simulation,  ASW  System  Evaluation  Tool 


19.  ABSIKACT  (Confimu  on  rtvtrst  if  ntctssary  and  idtniify  by  block  number) 

The  Antisubmarine  Warfare  System  Evaluation  Tool  (ASSET)  is  a  generic  high-level  antisubmarine  warfare  (ASW) 
modeling  tool,  designed  to  aid  ASW  personnel  in  ihc  development  and  refinement  of  ASW  top-level  warfare  requirements 
and  the  ASW  Master  Plan.  The  primary  objective  of  this  thesis  is  to  analyze  and  implement  the  improvements  suggested 
in  previous  evaluations  of  various  sub-areas  of  ASSET.  The  glimpse  raie  model  for  submarine  detection  used  in  ASSET 
hus  been  substituted  with  compound  Lambda-Sigma  jump  model.  There  is  a  different  target  radiated  frequency  in  each 
environmental  region.  Each  target  will  have  its  own  detection  rate  to  reflect  (he  differences  in  its  operating  characteristics. 
Multiple  engagements  between  platforms  are  used  to  eliminate  the  limitations  of  interaction  between  opponent  platforms. 
The  glimpse  rate  model  is  used  to  determine  dcicciion  opponuniiies  of  maritime  patrol  aircraft  (MPA)  and  (o  approximate  a 
continuous-looking  sensor  pattern,  A  different  criterion  of  selecting  search  probability  area  (SPA)  and  MPA  pairs  using  the 
ratio  of  MPA's  time  on-stalion  over  the  SPA  size  was  implemented.  The  feasibility  of  converting  current  ASSET  code  to 
CLOS  was  investigated.  In  addition,  part  of  the  code  was  converted  to  CLOS. 


p  tiij  I  bvti  PiVvi  I ; 


20.  DISTRffiunO.VAVAlLABlLnY  OF  ABSTRACT  T2I.  ABSTRACTS 

UNa-ASOTIED/UyLlMITED  Q]  SA-MEASRIT.  Q]  DDC USERS  UNCLASSIFIED 


22*.  NA.V1E  OF  REiPO.NSlDLE  INDIVIDUAL 

Yuh-Jeng  Lee 


cation 


22c.  OFFICE  SYMBOL 

CS/Le 


BD  FORM  1473,  84  MAR 


22b.  TELEPHONE  (Include  Aree 
Code) 

(408)  646-2361 


83  Al’K  cdiiiun  may  be  uicd  uniJ  exhauiied  SECURITY  CLASSIhCATION  OF TlflS  PACE 

AU  other  cJiiions  arc  obiolete  UNCLASSIFIED 


1 


Approved  for  public  release;  distribution  is  unlimited. 


Evaluation  and  Improvement 
of 

the  ASW  System  Evaluation  Tool 


by 

Peng-tso  Chang 

LT.  Commander,  Republic  of  China  Navy 
B.S..  Chinese  Naval  Academy,  1980 


Submitted  in  partial  fulfillment  of  the  requirements 
for  the  degree  of 

MASTER  OF  SCIENCE  IN  COMPUTER  SCIENCE 

from  the 


Author: 
Approved  by: 


NAVAL  POSTGRADUATE  SCHOOI. 
March  1992 


Department  of  Computer  Science 


ABSTRACT 


The  Antisubmarine  Warfare  System  Evaluation  Tool  (ASSET)  is  a  generic  high- 
level  antisubmarine  warfare  (ASW)  modeling  tool,  designed  to  aid  ASW 
personnel  in  the  development  and  refinement  of  ASW  top-level  warfare 
requirements  and  the  ASW  Master  Plan.  The  primary  objective  of  this  thesis  is 
to  analyze  and  implement  the  improvements  suggested  in  previous  evaluations  of 
various  sub-areas  of  ASSET.  The  glimpse  rate  model  for  submarine  detection 
used  in  ASSET  has  been  substituted  with  compound  Lambda-Sigma  jump  model. 
There  is  a  different  target  radiated  frequency  in  each  environmental  region. 
Each  target  will  have  its  own  detection  rate  to  reflect  the  differences  in  its 
operating  characteristics.  Multiple  engagements  between  platforms  are  used  to 
eliminate  the  limitations  of  interaction  between  opponent  platforms.  The  glimpse 
rate  model  is  used  to  determine  detection  opportunities  of  maritime  patrol 
aircraft  (MPA)  and  to  approximate  a  continuous-looking  sensor  pattern.  A 
different  criterion  of  selecting  search  probability  area  (SPA)  and  MPA  pairs 
using  the  ratio  of  MPA's  time  on-station  over  the  SPA  size  was  Implemented. 
The  feasibility  of  converting  current  ASSET  code  to  CLOS  was  investigated.  In 
addition,  part  of  the  code  was  converted  to  CLOS. 
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I.  INTRODUCTION 


A.  OVERVIEW  OF  ASSET 

The  Antisubmarine  Warfare  Systems  Evaluation  Tool  (ASSET)  is  a  generic 
high-level  ASW  modeling  and  analysis  tool.  It  is  designed  to  aid  ASW  personnel 
in  the  development  and  refinement  of  ASW  top-level  warfare  requirements  and 
the  ASW  Master  Plan.  Current  version  of  ASSET  (version  1.0)  was 
implemented  in  Macintosh  Allegro  Common  Lisp  (MACL)  version  1.3.2  and 
designed  for  use  on  a  Apple  Macintosh  II  computer. 

The  ASSET  is  used  to  model  open  ocean  scenarios  involving  submarines, 
maritime  patrol  aircraft  (MPA),  shore-based  command  and  data-fusion  centers, 
and  a  variety  of  passive  acoustic  sensors.  These  tactical  objects  interact  in  a 
panicular  geographic  location  against  an  analogous  enemy  force  structure.  The 
simulation  scenario  is  specified  by  a  user-supplied  architecture  that  determines 
all  aspects  of  environment,  command  control,  sensor  interaction,  and  platform 
maneuver.  This  structure  is  created  through  a  series  of  predefined  windows, 
each  dealing  with  a  specific  topic,  such  as  Geoplot  Edit  Window,  Track  Input 
Window,  and  Region  Build  Edit  Window. 

Once  the  parameters  are  determined  by  the  user,  the  scenario  is  repeated  as  a 
Monte  Carlo  simulation  to  produce  statistically  meaningful  measure  of 
effectiveness  (MOE).  Output  data  regarding  the  detection,  localization,  and 
prosecution  of  enemy  platforms  can  then  provide  a  quantitative  basis  for 
decisions  regarding  assessments  and  system  appraisals.  The  modular  nature  of 
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the  object-oriented  structure  of  ASSET  makes  expanding  the  scope  of  the 
simulation  possible. 

The  ASSET  is  programmed  in  Macintosh  Allegro  Common  Lisp,  an  object- 
oriented  programming  language,  which  allows  related  elements  in  the  simulation 
to  inherit  data  and  functionality  via  hierarchical  class  structure.  The  simulated 
objects  in  ASSET  can  be  broken  into  three  groups;  (1)  command,  control,  and 
communication  (C3)  objects;  (2)  acoustic  radiators;  and  (3)  ASW  detectors.  The 
C3  objects  consist  of  level  1  commands,  level  2  commands,  fusion  centers, 
antisubmarine  warfare  operational  centers  (ASWOCs),  Submarine  Operating 
Authorities,  and  communication  satellites.  A  more  detailed  description  of 
ASSET  and  its  capabilities  can  be  found  in  the  ASSET  Technical  Documentation 
and  User's  Manual  [R1  90). 

Weaknesses  of  the  ASSET  that  have  been  identified  include  the  following: 
unable  to  use  lower-level  models,  only  part  of  ASW  problem  is  considered, 
difficult  for  the  novice  to  use,  excessively  simple  acoustic  and  maritime  patrol 
aircraft  (MPA)  modeling,  and  slow  execution  time. 

B.  PREVIOUS  WORK  AT  NPS 

Several  evaluations  on  the  various  sub-models  of  ASSET  have  been 
conducted  at  the  Naval  Postgraduate  School.  [CA  91]  examined  and  evaluated  the 
acoustic  detection  modeling  in  ASSET.  An  alternate  detection  model  has  been 
suggested  that  incorporates  aspects  of  the  glimpse  rate  model  and  the  more 
standard  Lambda-Sigma  (A,-o)  model. 

(SH  91]  evaluated  the  MPA  detection  and  allocation  models  utilized  by 
ASSET.  It  also  proposed  a  potential  long  term  course  of  action  for  evolving 
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ASSET  into  a  comprehensive,  flexible,  and  simple-to-use  tool  for  top  level  ASW 
appraisals  and  assessments. 

[VE  91]  investigated  the  mathematical  development  of  the  tracking 
algorithm.  Some  possible  modifications  to  the  existing  Kalman  filter  tracking 
algorithm  were  suggested  to  better  match  the  data  input  and  increase 
computational  efficiency. 

C.  OBJECTIVES 

The  primary  objective  of  this  thesis  is  to  implement  and  verify  the 
improvements  suggested  in  previous  NFS  evaluations  of  various  sub-areas  of 
ASSET.  We  want  to  test  whether  the  suggestions  are  feasible.  A  second 
objective  is  to  improve  the  overall  efficiency  of  the  ASSET  run-time  system. 
The  feasibility  of  converting  current  Object  Lisp  code  to  Common  Lisp  Object 
System  (CLOS)  is  also  being  investigated. 

D.  ORGANIZATION  OF  THESIS 

Chapter  II  provides  an  overview  of  object-oriented  programming  language 
(OOPL)  in  simulation  application.  W'e  review  the  basic  characteristics  and 
benefits  of  OOPL  and  some  applications  in  the  object-oiiented  simulation  field. 
Chapter  III  describes  the  current  detection  modeling  and  possible  improvements, 
followed  by  a  description  of  implementation  changes.  Chapter  IV  contains  the 
description  of  the  current  MPA  modeling,  possible  improvements,  and  the  details 
of  implementation.  Chapter  V  presents  the  work  of  converting  current  version 
to  CLOS.  Finally,  Chapter  VI  is  the  conclusion  and  summary  along  with  a  brief 
statement  of  recommendations  for  future  work.  Appendix  A  contains  the  source 
code  of  rei'ised  and  new  modules  of  ASSET. 
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II.  OBJECT-ORIENTED  SIMULATION 


In  discrete  event  simulations,  we  normally  simulate  the  behavior  of  various 
objects  that  may  be  mathematical  or  statistical  constructs,  or  in  most  cases, 
physical  and  observable  entities.  For  example,  in  a  manufacturing  factory,  the 
objects  are  the  machines  of  interest,  the  types  of  parts  that  are  to  be  produced, 
and  the  operations  needed  to  complete  their  manufacture.  Such  objects  are 
usually  the  central  focus  of  the  simulation  studies  and  are  easy  to  identify.  The 
purposes  of  simulations  are  to  find  a  convenient  means  of  modeling  the  objects 
and  to  manipulate  and  control  their  behaviors,  including  the  interactions  between 
the  objects. 

In  a  simulation,  the  user  creates  computer  models  of  the  various  elements  of 
the  simulation,  describes  how  they  will  interact  with  one  another,  and  sets  them 
moving.  This  is  similar  to  typical  object-oriented  programming,  in  which  the 
user  creates  the  entities  in  the  universe  for  the  program,  describes  how  they  will 
interact  with  one  another,  and  finally  sets  them  in  motion.  Due  to  the 
similarities,  there  are  many  advantages  in  using  the  object-oriented  paradigm  for 
discrete  event  simulations. 

A.  BASIC  CONCEPTS  OF  OBJECT-ORIENTED  COMPUTING 
1 .  Objects,  Classes,  and  Inheritance 

Object-oriented  languages  combine  the  descriptions  of  data  and 
procedures  within  a  single  entity  called  an  object.  An  object  consists  of  well 
defined  set  of  variables  (called  slots  within  the  object)  that  can  be  manipulated 
only  by  a  set  of  methods  defined  exclusively  for  that  purpose.  A  class  is  a 
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description  of  one  or  more  similar  objects.  It  is  used  to  specify  the  implemented 
data  structures  and  the  operations  that  can  be  performed. 

The  variables  making  up  an  object  can  be  divided  into  two  kinds:  class 
variables  and  instance  variables.  A  class  variable  is  the  one  we  associate  to  the 
class  itself  and  is  shared  in  both  name  and  value  by  all  instances  of  a  class.  An 
instance  variable  is  shared  in  name  only  by  all  instances  of  a  class. 

A  method  is  a  procedure  or  operation  that  is  defined  for  an  object.  A 
message  is  the  instruction  sent  to  an  object  to  perform  one  of  its  methods. 
Ideally,  the  only  way  to  access  any  of  the  variables  making  up  an  object  is  by 
sending  the  object  a  message.  In  this  way,  an  object  is  said  to  be  encapsulated  in 
that  its  internal  structure  may  be  modified  without  affecting  us ’r- written  code 
that  accesses  the  object. 

Inheritance  can  be  defined  simply  as  code  sharing  mechanism.  It  allows 
a  new  class  to  be  defined  based  upon  the  definition  of  an  existing  class  without 
having  to  copy  all  the  existing  code.  A  child  class  (or  subclass)  inherits  all  the 
variables  and  methods  defined  for  its  parent  class  (or  superclass).  Inheritance  is 
usually  transitive,  so  a  class  can  inherit  features  from  superclasses  many  levels 
above.  Exceptions  can  be  handled  at  any  level  by  altering  the  values  of  the 
affected  variables. 

Inheritance  means  that  the  behavior  and  data  as-ociated  with  child 
classes  are  always  extensions  of  ilie  properties  associated  with  parent  classes.  A 
subclass  must  have  all  the  properties  of  the  parent  class,  and  others  as  well.  On 
the  other  hand,  since  a  child  class  is  a  more  specialized  (or  restricted)  form  of 
the  parent  class,  it  is  also,  in  a  certain  sense,  a  contraction  of  the  parent  type. 
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2.  Encapsulation 

Data  abstraction  is  a  methodological  approach  to  problem  solving  in 
programming  where  information  is  consciously  isolated  and  hidden  in  only  one 
part  of  a  program.  In  particular,  the  programmer  develops  a  series  of  abstract 
data  types.  Each  abstract  data  type  can  be  viewed  as  having  two  faces.  From  the 
outside,  a  client  (user)  of  an  abstract  data  type  sees  only  a  collection  of 
operations  that  together  characterize  the  behavior  of  the  abstraction.  On  the 
other  side  of  the  interface,  the  programmer  defining  the  abstraction  sees  the  data 
variables  that  are  used  to  maintain  the  internal  state  of  the  object.  For  example, 
in  a  graphics  package,  a  graphical  object  may  be  defined  to  have  interval  values 
of  position,  shape,  and  extent;  on  the  other  hand,  operations  such  as  moving  and 
combining  objects  are  all  that  the  user  needs  to  change  the  internal  values  of  an 
object.  Once  declared,  the  graphical  object  then  becomes  a  data  type  within  the 
language.  In  a  data  base  application,  a  record  becomes  a  general  data  type  and  its 
operations  may  include  storage  and  retrieval  of  the  records. 

3.  Polymorphism 

Polymorphism  can  be  defined  as  allowing  a  method  to  operate  on 
objects  of  different  data  types.  True  polymorphism  allows  the  various 
operations  to  be  created  and  tested  independently  of  one  another.  Adding  a  new 
version  of  the  operation  has  no  effect  upon  the  existing  code  at  all. 
Polymorphism  allows  the  definition  of  flexible  software  elements  amenable  to 
extension  and  reuse. 

.Simple  polymorphism  allows  each  of  different  classes  to  have  its  own 
implementation  of  an  operation.  Multiple  polymorphism  allows  each  class  to 


6 


have  several  operations  with  the  same  name.  The  proper  operation  is  chosen 
based  upon  the  arguments  provided. 

4.  Dynamic  Binding 

Many  benefits  of  object-oriented  programming  result  from  its  dynamic 
binding.  When  you  call  a  function  in  a  language  such  as  C,  the  compiler  and  the 
linker  cooperate  to  generate  a  call  to  a  physical  address.  While  this  is  very 
efficient,  one  must  take  care  to  associate  the  function  with  the  appropriate  data 
structures.  Strongly  typed  languages  attempt  to  catch  mismatched  data  types  at 
compile  time.  Others  do  a  poorer  job  of  catching  these  errors,  and  sometimes 
the  problems  are  not  detected  until  the  output  starts  looking  strange. 

In  object-oriented  programming  languages,  the  programme*’  is  relieved 
of  the  burden  of  calling  the  right  method  with  the  specific  data  structure. 
Instead,  the  programmer  uses  a  generic  name  for  the  function,  and  the  receiving 
object  looks  up  for  the  proper  method.  This  run-time  binding  is  sometimes 
called  "late  binding"  and  has  several  advantages.  In  the  programming  world,  late 
binding  means  that  references  are  symbolic  and  methods  can  be  compiled  without 
re-compiling  all  its  callers.  The  same  symbolic  names  are  used,  despite  the  type 
of  object.  Finally,  a  single  message  can  invoke  several  methods.  This  is  known 
as  polymorphic  behavior  and  it  allows  code  to  be  written  independent  of  the 
receiver. 

For  simulation,  late  binding  means  that  the  specific  machine  needed  for 
the  part  can  be  determined  when  the  part  finishes  its  previous  operation.  It 
doesn't  have  to  be  established  at  the  start  of  the  simulation.  The  time  to  complete 
the  part  can  also  depend  on  the  current  set  of  resources  and  the  state  of  the 
system.  For  example,  consider  a  conveyor  handeling  packages.  How  those 
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packages  occupy  space  on  the  conveyor  is  determined  at  the  time  they  are  placed 
there,  not  in  some  predetermined  set  of  fixed-size  bins. 

B.  MODELING  PERSPECTIVES 

Modeling  with  objects  focuses  on  what  the  objects  are  in  the  real  system, 
instead  of  what  the  objects  in  a  simulation  language  represent.  In  a  simulated 
system  the  general  entities  can  be  grouped  as  classes,  and  specific  entities  should 
be  objects.  If  the  objects  appear  to  share  some  common  characteristics,  then  they 
can  be  subclass  of  a  more  general  class.  Methods  should  be  defined  for  specific 
operations  if  needed.  Division  of  objects  into  classes,  recognition  of  methods, 
and  the  organization  of  hierarchies  form  the  basic  approach  to  object-oriented 
modeling. 

1.  Potential  Benefits  of  Object-oriented  Approach 

a.  Modular  Implementation 

Object-oriented  programming  enables  the  programmer  to  define  an 
organization  of  classes  that  models  the  relationships  among  the  various  kinds  of 
objects.  The  programer  can  define  classes  that  serve  as  building  blocks.  Each 
individual  aspect  of  structure  and  behavior  is  abstracted  and  defined  separately. 
The  programmer  then  creates  new  classes  that  inherit  the  desired  combination  of 
building  blocks.  Therefore,  object-oriented  approach  allows  for  a  modular 
implementation. 

b.  Conveniently  Extensible 

The  Object-oriented  simulation  is  conveniently  extensible.  Using  a 
set  of  classes  with  documented  structure  and  behavior  as  building  blocks,  users 
can  create  new  classes  and  add  customized  behavior.  Through  function  and 
operator  overloading,  old  symbols  take  on  additional  meaning.  Inheritance 
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permits  new  objects  to  be  defined  from  existing  ones,  only  the  differences  need 
to  be  noted.  Old  models  are  reusable  because  their  methods  and  objects  continue 
to  be  useful. 

c.  Reduced  Code  Size 

In  side-by-side  comparisons  code  written  in  object-oriented 
programming  languages  arc  substantially  smaller  than  that  written  in  procedural 
languages  [CO  87],  The  reduced  code  size  means  that  a  single  person  can  manage 
more  tasks.  In  the  simulation  of  complex  systems,  building  larger  and  more 
realistic  models  are  possible  without  an  increase  in  manpower. 

d.  Natural  Representation  of  Objects 

Objects  in  most  simulations  tend  to  be  physical  and  real,  and  can 
easily  be  represented  pictorially.  Therefore,  object-oriented  simulation  models 
often  have  a  natural  pictorial  (iconic)  representation  and  are  easily  animated. 
The  user  can  often  translate  directly  his  simulation  model  into  an  animated 
simulation  without  additional  conceptual  changes. 

Since  the  objects  contain  their  own  functionality,  intelligence  can  be 
built  directly  into  this  functionality  using  the  machinery  of  artificial  intelligence 
and  expert  systems.  In  addition,  objects  provide  a  natural  basis  for  concurrency. 
The  idea  is  that  each  object  could  be  assigned  to  its  own  processor  and  run 
independently  until  some  form  of  coordination  is  needed.  Although  it  isn't  clear 
exactly  what  form  the  coordination  should  take,  there  is  a  natural  division  among 
the  simulation  components  when  viewed  as  objects. 

e.  Reusable  Software  Systems 

By  reducing  the  interdependency  among  software  components, 
object-oriented  programming  permits  the  development  of  reusable  software 
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systems.  Such  components  can  be  created  and  tested  as  independent  units, 
isolated  from  other  software  components. 

2.  Disadvantages 

a.  Run-time  Cost 

Dynamic  binding  is  flexible  but  with  some  run-time  cost.  Dynamic 
method  binding  requires  the  execution  of  some  run-time  mechanism  to  match  a 
method  to  a  message.  The  cost  of  dynamic  typing  involves  a  run-time  search  to 
discover  the  code  to  execute  every  time  operation  is  used  on  a  data  value. 

b.  Memory 

Object-oriented  environments  (like  CLOS)  require  machines  with 
lots  of  RAM.  In  addition,  the  use  of  any  software  library  frequently  imposes  a 
size  penalty  over  the  use  of  systems  especially  constructed  for  a  specific  project. 
With  the  advance  in  hardware  technology,  however,  this  has  become  less  serious. 

c.  Program  Complexity 

Although  object-oriented  programming  is  often  touted  as  a  solution 
to  the  problem  of  software  complexity,  overuse  of  inheritance  often  replaces  one 
form  of  complexity  with  another.  Understanding  the  control  flow  of  a  program 
that  uses  inheritance  may  require  several  multiple  scans  up  and  down  the 
inheritance  graph. 

Some  object-oriented  languages  require  extensive  class  libraries  be 
understood  before  becoming  proficient.  This  increases  the  learning  time  and 
forces  users  to  become  more  dependent  on  documentation  and  high-level 
debugging  tools. 
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d.  No  Standard  Model 

Different  object-oriented  approaches  use  different  terminology  to 
define  similar  concepts.  There  is  no  single,  agreed  upon  standard  for  object- 
oriented  model. 

C.  SIMULATIONS  WITH  OBJECT-ORIENTED 
IMPLEMENTATIONS 

Many  object-oriented  systems  are  suitable  for  performing  simulations.  None 
have  clearly  been  proven  to  be  generally  superior  as  research  into  understanding 
object-oriented  programming  continues.  Among  the  widely  used  languages, 
Simula  was  designed  as  a  simulation  language,  Smalltalk  (GO  83]  contains  an 
extensive  set  of  classes  and  methods  to  support  simulation,  and  C-n-  [ST  87]  was 
designed  to  be  applied  to  simulation.  Uses  of  Smalltalk  as  a  simulation  language 
and  environment  can  be  found  in  (  GO  83],  and  in  the  tutorial  of  (KN  86,  KN 

87] . 

There  are  several  Lisp- based  object-oriented  languages  with  mechanisms  for 
abstraction,  polymorphic  typing,  and  dynamic  binding.  These  include  Flavors 
(which  supports  a  more  primitive  function  than  an  object),  CommonLoops,  and 
New  Flavors.  A  Lisp  based  object-oriented  simulation  system  is  described  in  [ST 

88] ,  DEVS-Scheme  is  an  implementation  of  DEVS  for  hierarchical,  modular 
system  within  an  object-oriented  framework  (K1  87].  DEVS  is  a  simulation 
formalism  developed  in  [ZE  84). 

There  are  also  several  hybrid  object-oriented  s  - stems  that  combine  the 
object-oriented  approach  with  traditional  procedural  features.  For  example, 
Objective-C  [CO  87]  adds  objects,  similar  to  Smalltalk,  to  the  definition  of  C. 
With  the  similarity  between  Smalltalk  and  Objective-C,  a  translator  has  been 
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recently  introduced  to  convert  Smalltalk  into  Objective-C  [SC  87].  Thus 
Smalltalk  could  be  used  for  design  and  prototyping,  and  later  converted  to  C 
through  Object-C  for  efficient  execution.  Actor  [DU  86]  provides  a  Smalltalk 
type  environment  with  Pascal-like  procedural  programming  and  artificial 
intelligence  features. 

Object-oriented  systems  provide  a  practical  approach  for  those  who  design 
simulation  software.  Object-oriented  languages  provide  a  natural  framework  for 
development.  The  information  hiding  and  abstraction  facilities  make  it  easy  to 
develop  and  maintain  complex  software  components.  The  extensible  platform  is 
an  attractive  way  to  add  new  concepts  and  features  to  an  existing  simulation.  It  is 
just  a  matter  of  time  before  existing  simulation  languages  attempt  to  exploit 
various  aspects  of  the  powerful  features  of  the  object-oriented  systems. 
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ni.  ACOUSTIC  DETECTION  MODELING 


A  detection-rate  model  is  used  in  ASSET  version  1.0  to  model  acoustic 
systems  in  which  detection  opportunities  occur  over  continuous  intervals  of  time. 
It  is  used  for  modeling  the  detection  capabilities  of  submarines,  maritime  patrol 
aircraft  (MPA),  fixed  area  sensor  (FAS),  mine  fields,  and  trip  wires. 

The  first  section  gives  a  brief  description  of  die  existing  ASSET  model  and  its 
possible  problems.  The  second  sectia  reviews  the  alternate  model  suggested  in 
[CA  91  ]  and  provides  some  addition:'.’  ideas  on  improving  the  current  model.  The 
implementation  details  for  improvement  are  discussed  in  the  third  section. 

A.  EXISTING  ASSET  MODEL 
1.  Glimpse  Rate  Model 

Because  .u'  computational  limitations,  ASSET  uses  a  discrete  glimpse 
rate  to  model  :he  c.  .,.tinu:.us  process  of  passive  acoustic  detection.  By  drawing  an 
exponentially  disiributed  random  variable  for  the  time  of  the  next  glimpse  (i.e., 
detection  opportunity),  run  time  for  complex  simulations  can  be  held  to  an 
acceptable  level. 

In  the  glimpse  rate  model,  detection  opportunities  occur  at  times 
determined  by  a  Poisson  process  with  user  specified  rate  g.  When  an  opportunity 
occurs,  the  target  is  detect>yJ  if  and  only  if  signal  excess  (SE(r))  is  greater  than 
zero.  The  details  of  glimpse  rate  ti.odel  used  in  ASSET  can  be  found  in  (RI 90]. 

The  glimpse  rate  is  a  user  input  which  greatly  affect  the  results  of  a  run. 
(CA  91]  points  out  that  there  is  no  good  guidance  for  selecting  the  glimpse  rate  g. 
If  g  is  small,  the  time  between  glimpses  can  be  quite  long.  Rapid  changes  in  the 
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mean  signal  excess  may  be  missed,  causing  searchers  to  have  little  to  no  success 
detecting  other  platforms.  If  g  is  too  large,  the  time  between  glimpse  can  become 
small  enough  that  the  model  approaches  continuous  glimpsing,  resulting  in  a 
unrealistically  high  cumulative  detection  probability  CDP(/).  This  happens 
because  ASSET  draws  an  independent  fluctuation  value  for  each  look,  possibly 
underestimating  the  expected  correlation  between  signal  excess  values  at  closely 
spaced  times. 

Another  major  drawback  is  that  currently  there  is  only  one  radiated 
frequency  for  each  acoustic  platform  in  ASSET.  In  addition,  there  is  a  single 
mean  glimpse  interval  by  all  searchers  regardless  of  the  searcher  type  or  mission. 

2.  IiUvraction  between  Acoustic  Platforms 

During  the  period  of  a  hostile  encounter,  a  platform  can  engage  with  (or 
be  engaged  by)  <.  nly  one  opposing  platform.  If  a  platform  is  engaging  with  an 
opponent,  then  it  will  not  be  detectable  by  any  other  platform. 

B.  SUGGESTED  IMPROVEMENTS  TO  ASSET 
1.  Compound  Lambda-Sigma  Jump  Process 

The  Lambda-Sigma  jump  process  as  described  in  [AVA  91]  is  a  model 
where  the  times  of  random  fluctuations  in  signal  excess  occur  according  to  a 
Poisson  process  with  rate  X  per  unit  time.  The  amount  of  the  fluctuation  at  time  t 
(D(/))  is  drawn  independently  from  a  normal  distribution  with  mean  zero  and 
standard  deviation  a.  The  deviation  value  remains  constant  until  the  next 
fluctuation  occurs  in  the  Poisson  process,  at  which  time  a  new  deviation  is  drawn, 
independent  of  previous  values.  The  deviation  step  function  is  added  to  the  mean 
signal  excess  curve,  resulting  in  a  discontinuous  signal  excess  curve.  Like  the 
glimpse  rate  model,  a  detection  occurs  if  signal  excess  is  nonnegative. 
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The  advantages  of  using  the  Lambda-Sigma  jump  process  in  ASW 
analyses  are: 

(1)  It  models  random  fluctuations  of  signal  strength  in  a  handy  and  adaptable 

manner. 

(2)  It  captures  the  idea  of  correlation  in  signal  excess  between  closely  spaced 

times. 

(3)  For  simple  geometries,  close-form  expressions  for  CDP  (t)  are  available. 

(4)  It  is  extremely  easy  to  simulate. 

a.  Global  Lambda-Sigma  Jump  Process 

The  ocean  fluctuation  value  determined  by  a  Lambda-Sigma 
stochastic  process  is  used  to  simulate  global  area  environmental  correlation  for  all 
targets.  This  value  is  used  as  one  of  the  variables  in  computing  the  signal  excess. 

b.  Indi vidual  Lam  bda-Sigma  Jump  Process 

To  reflect  the  acoustic  environmental  differences  among  the  various 
searchers'  local  environment  regions,  an  individual  Lambda-Sigma  value  is  also 
used  as  one  of  the  variables  in  computing  the  signal  excess. 

c.  Total  Signal  Excess 

Total  signal  excess  is  the  sum  of  the  mean  signal  excess  plus  the 
global  and  individual  Lambda-Sigma  fluctuation  values. 

2.  Target's  Radiated  Frequency 

The  most  detectable  radiated  frequency  of  a  target  will  be  selected  from  a 
list  of  environment-frequency  pairs  according  to  the  environmental  region  it  is  in. 
The  frequency  is  then  used  to  determine  propagation  loss  and  sweep  width  in 
computing  the  signal  excess.  This  will  allow  a  different  target  radiated  frequency 
in  each  environmental  region. 
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3.  Target's  Detection  Rate 

Each  target  will  have  its  own  detection  rate  to  reflect  the  differences  in  its 
operating  characteristics.  Consequently,  each  searcher-target  pair  will  have  its 
own  mean  time  (between  glimpse)  for  scheduling  detection  opportunity  events. 

4.  Multiple  Engagements  between  Submarines 

When  a  platform  is  engaging  with  an  opponent,  it  is  still  detectable  to 
other  opponent  searchers,  but  it  will  not  detect  any  more  platforms  before  it 
finishes  the  engagement.  Thus  it  is  possible  for  a  platform  to  be  delected  and 
attacked  while  it  is  engaging  a  third  platform.  To  simplify  the  calcuiation  of  the 
probability  of  a  kill  in  this  complex  engagement,  the  events  of  engagement  are 
sequentially  processed.  The  probability  of  a  kill  is  treated  independently  among 
all  searcher-target  pairs  involved  in  this  multiple  engagement.  The  system  will 
check  that  either  the  searcher  or  target  is  active  before  processing  engagement 
events  and  platform  kill  events. 

5.  Responding  False  Alarm  Target 

When  a  false  target  appears  during  search,  the  searcher  will  send  a 
detection  report  and  will  continue  the  search  operations.  Also,  the  platform  having 
a  false  detection  will  still  be  detectable  to  other  unengaged  platforms. 

C.  IMPLEMENTATION  OF  IMPROVEMENT  IDEAS 

In  this  section  we  describe  implementation  ideas  of  the  improvements 
mentioned  above.  The  changes  can  be  classified  into  three  categories:  data 
structures,  modified  functions,  and  new  functions.  They  are  categorized  by  their 
physical  locations  in  the  application  package  distributed  by  Metron. 
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1.  Structure  File 

a.  PlatformacousticParms 

The  structure  Platformacousticparms  specifies  acoustic  parameters 
of  a  platform.  The  following  data  fields  are  modified: 

individualSigma  Sigma  for  individual  Lambda-Sigma 

process 

LmabdaSig.TiaValue  Individual  fluctuation  of  target's  signal 

excess 

meaninterval  Mean  time  between  fluctuation 

detectableincerval  Target's  detectable  interval  (1/Lambda  ) 

env-f  req-list  List  of  target's  ladiated  frequencies  in  each 

region 

2.  RB-Instailer  File 

a.  New  Radio  Buttons 

Two  new  radio  buttons  are  added  for  input  of  the  compound 
Lambda-sigma  jump  process: 

lambdaSigtnaProcess-  Radio  button  for  input  of  Lambda-Sigma 

rb  parameters 

env-freq-rb  Radio  button  for  input  of  environment- 

frequency  pairs  of  a  acoustic  platform 

3.  Class  Umpire 

The  class  Umpire  describes  the  variables  and  methods  to  peiform  the 
setup  of  platform  detection  and  engagement  event  schedules.  It  is  modified  as 
follows: 


17 


a  Data  Field  Deleted 


FOMglimpse Interval 

b.  New  Data  Fields 

global  Sigma  Sigma  for  global  Lambda-Sigma 

process 

globalLambdaSigmaValue  Global  fluctuation  for  signal  excess 
globalMeanInterval  Mean  time  between  global  fluctuation 

c.  New  Method 

updateGlobalLambdaSiginaValue  :  Updates  global  Lambda- 
Sigma  value  and  schedules  for  next  update. 

d.  Modified  Functions 

setupForReplicatioir :  Initializes  globalLambdaSigmaValue  and 
schedules  for  next  update. 

aatNextDetectionOpportunity :  Allows  engaged  target  being 
detectable  to  searchers  that  are  not  yet  engaged.  It  also  checks  range 
when  scheduling  next  detection  opportunity  for  engaged  platforms. 
processOetect  ionOpportur4ix;y :  Engaged  platforms  will  not 
detect  any  other  target. 

setNextDetectionOpportunity :  Allows  engaged  target  being 
detectable  to  searchers  that  are  not  yet  engaged.  It  also  checks  rarge 
when  scheduling  next  detection  opportunity  for  engaged  platforms. 
processDetect  ionOpportunity :  Engaged  platforms  will  not 
detect  any  other  target.  Engaged  platforms  is  detectable  to  unengaged 
searchers. 
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gatDetectionOppRdsult :  Determines  target's  most  detectable 
radiated  frequency.  Uses  global  and  individual  Lambda-Sigma  value  to 
adjust  mean  signal  excess. 

processPlat£ormKill :  Checks  if  the  platform  had  been  killed 
already. 

aatNextFalaeAlarm :  Checks  if  the  platform  is  aedve  before  it  sets 
up  another  false  alarm. 

procasaFalseAlarm:  Checks  if  the  platform  is  active  before  it 
processes  the  false  alarm. 

4.  ClEiS  EnvironmentManager 

The  class  EnvironmentManager  contains  the  acoustic  characteristics 
information  for  each  environmental  region,  and  defines  the  methods  of  getting 
propagation  loss  and  sweep  width.  Its  modifications  are  described  as  follows: 

a.  Data  Fields  Deleted 
sigma 

b.  Modified  Functions 

g«tPL-AN-Sign\c. :  Target's  most  detectable  frequency  is  used  to 
determine  propagation  loss  and  ambient  noise. 

getSweepWidth :  Target’s  most  detectable  frequency  is  used  to 
determine  sweep  width. 

5.  Class  AcousticPIatforms 

The  code  that  describes  the  variables  and  operations  of  an  acoustic 
platform  is  changed  as  follows; 

a.  Data  Fields  Deleted 

glimpse  Interval  Now  using  target's  detectable  interval 
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b.  New  Data  Fields 

individualSigma  Sigma  of  individual  Lambda-Sigma 

process 

iambdaSigmaValue  Individual  fluctuation  for  signal  excess 
meaninterval  Mean  time  between  updating  individual 

Fluctuation 

detectableinterval  target's  detectable  interval  used  to  schedule 

next  glimpse 

c.  Modified  Functions 

editEnv-Freq :  Inputs  environment-frequency  pairs  of  a  platform. 
setForStart :  Initializes  env-freq-list,  Lambda-Sigma  process 
parameters,  target's  detectable  interval  and  schedules  event  for  next 
update  of  individual  Lambda-Sigma  value. 

satForStart :  Initializes  env-frcq-list,  Lambda-Sigma  process 
parameters,  target's  detectable  interval  and  schedules  event  for  next 
update  of  individual  Lambda-Sigma  value. 

d.  New  Functions 

updateLambdaSigmaValue :  Updates  Lambda-Sigma  value  and 
schedule  event  for  next  update  of  individual  Lambda-Sigma  value. 

6.  Changes  on  User  Inputs 

'fhis  section  explains  the  changes  on  some  necessary  inputs  that  the  user 
needs  to  specify  in  setting  up  the  simulation.  The  details  of  setting  up  a  whole 
simulation  run  can  be  found  in  [RI  90).  The  user  now  should  enter  Global  Mean 
Interval  (1  A)  and  Global  Sigma  during  setting  up  Umpire  Parameters.  This  can  be 
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done  by  choosing  the  Edit  Umpire  Parameters  menu  item  under  Simulation  menu, 
which  displays  the  Simulation  Umpire  Edit  window  shown  in  Figure  1. 


Figure  1 .  Simulation  Umpire  Edit  Windov 


The  Edit  Submarine  window  in  Figure  2  contains  added  radio  buttons: 
Lambda-Sigma  Process,  Detectable  Interval,  and  Frequency  Environment  Pairs. 
To  create  or  edit  a  submarine  object,  the  user  needs  to  specify  appropriate 
parameters  in  Lambda-Sigma  Process  Edit  window  shown  in  Figure  3.  For  each 
environmental  region  the  user  must  enter  a  most  detectable  frequency  of  the 
platform.  The  Frequency  Environment  Edit  window  is  shown  in  Figure  4. 
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Figure  2.  Edit  Submarine  Window 
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IV.  MPA  MODELING 


In  this  chapter  we  review  ASSETs  current  MPA  detection  rate  model  and 
allocation  model,  including  discussion  of  the  disadvantages  and  suggested 
improvements.  This  is  followed  by  a  description  of  the  implementation  details. 

A.  CURRENT  MPA  DETECTION  RATE  MODEL 

The  MPA  acoustic  detection  model  uses  a  detection  rate  scheme  to  decide 
whether  a  detection  is  made.  This  model  is  derived  from  the  passive  sonar 
equation  and  the  random  search  formula. 

'Fhe  MPA  is  assumed  to  lay  buoys  in  a  pattern  to  provide  uniform  coverage  of 
the  designated  search  region.  When  an  MPA  is  cued  to  a  search  region,  the  search 
area  is  the  86%  containment  region  provided  by  the  tracker-correlator.  If  the  MPA 
is  not  cued,  then  the  search  region  has  user-input  size,  and  is  randomly  located  in 
the  user-input  area  search  region.  If  a  target  is  in  the  designated  MPA  search 
region  at  the  beginning  of  the  MPA's  time  on-station,  then  it  is  assumed  that  the 
target  remains  in  the  search  region  for  the  entire  MPA's  time  on-station.  And,  if 
the  target  is  outside  the  search  region  at  the  beginning  of  the  MPA's  time  on- 
station,  then  it  is  assumed  that  the  target  remains  outside  the  search  region  for  the 
entire  MPA's  time  on-station. 

1.  Detection  Rate  Calculation 

The  ASSET  calculates  a  detection  rate  for  each  target  that  lies  within  the 
search  region  when  the  MPA  arrives  on-station.  The  detection  rate  is  the  ratio  of 
the  area  searched  per  unit  time  over  the  total  area  of  the  search  region.  ASSET 
uses  a  constant  detection  rate  7  given  by: 
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y  =  NVW/As 


where  N  is  the  number  of  sonobuoy  channels  processed,  V  is  the  average  target 
speed, W  is  the  acoustic  sweep  width  of  a  single  sonobuoy,  and  As  is  the  search 
area. 

The  number  of  sonobuoy  channels  processed,  N,  is  the  smaller  of 
following  two  values:  user-selected  number  of  buoys  per  search,  or  the  user- 
selected  number  of  buoys  that  can  be  processed. 

The  single  buoy  acoustic  sweep  width,  W.  is  calculated  as  twice  the 
maximum  detection  range,  Rmax,  the  maximum  range  at  which  the  adjusted  or 
actual  figure  of  merit  (FOM)  is  equal  to  the  propagation  loss  from  the  user-entered 

Proploss  Table.  Mean  FOM  is  given  by: 

FOM  =  SL  -  NL  +  D1  -  DT. 

where  SL  is  the  target  radiated  source  level,  NL  is  the  total  noise  level  (self-noise 
+  ambient  noise),  DI  is  directivity  index  of  the  receiver,  DT  is  the  detection 
threshold  or  recognition  differential. 

Actual  FOM  is  obtained  by  adding  a  value  representing  the 
environmental  uncertainty  correction  to  the  mean  FOM.  The  uncertainty 
correction  is  a  normally  distributed  random  variable  with  a  mean  of  zero  and  a 
standard  deviation  of  9  dB.  ASSET  generates  a  single  uncertainty  correction,  and 
thus  a  single  Rmax.  that  is  used  for  the  entire  MPA  search  period.  SL,  NL,  DI, 
and  DT  are  user-entered  parameters  for  the  particular  MPA-submarinc  pair. 

2.  Determination  of  Detection 

After  the  detection  rates  are  calculated  for  each  submarine  that  is  within 
the  search  region  when  the  MPA  arrives,  ASSET  sums  each  of  these  rates  with  the 
user-entered  false  alarm  rate  to  obtain  a  collective  contact  rate  (F).  The 
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probability  distribution  for  the  time  to  initial  contact  is  an  exponential  distribution 
with  a  rate  equal  to  the  collective  contact  rate.  ASSET  generates  an  exponential 
random  number  t  that  represents  the  time  when  the  MPA  detects  a  real  or  false 
target.  The  exponential  random  number  is  obtained  from 

t  =  -ln(U[0.1J)/r 

where  U(0,1]  is  a  unifonn  random  number  between  zero  and  one.  A  detection  is 
reported  at  time  t  if  t  is  less  than  or  equal  to  the  total  time  that  the  MPA  is  on 
station. 

To  determine  if  a  target  in  the  search  region  is  detected,  ASSET  stacks 
the  detection  and  false  alarm  rates  end  to  end.  ASSET  then  draws  another  uniform 
random  number,  on  the  interval  (0,  Ht  to  determine  whether  a  submarine  (or  false 
alarm)  will  be  reported  as  the  contact.  ASSET  only  allows  one  detection  (real  or 
false)  per  MPA  mission. 

3.  Evaluations  of  Current  Model 

(SH  91]  has  indicated  several  drawbacks  in  current  ASSETs  MPA 
model.  Targets  outside  of  the  search  region  at  the  beginning  of  MPA’s  on  station 
are  excluded  the  search.  So  an  MPA  can  not  detect  targets  crossing  the  search 
region  boundary  during  its  searching. 

The  ability  of  an  MPA  to  detect  and  kill  a  target  is  affected  by  the  MPA 
and  submarine  false  alarm  rate.  Apparently  due  to  the  limitation  that  engaged 
platforms  are  not  detectable  by  opponent  platforms,  even  extremely  small 
submarine  false  alarm  rates  significantly  reduce  the  detectability  of  the  submarine. 

In  addition,  currently  only  a  single  contact  is  performed  per  MPA 
mission,  even  if  that  contact  is  a  false  alarm  and  the  MPA  has  sufficient  time  to 
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remain  on-station.  After  the  MPA  completes  one  engagement,  no  further  search 
is  conducted. 

B.  ALLOCATION  MODEL 

In  ASSET  each  antisubmarine  warfare  operational  center  (ASWOC)  has  a 
specific  non-overlapping  ocean  area  of  responsibility  and  several  assigned  MPA 
squadrons,  both  defined  by  the  user.  Based  upon  fused  submarine  picture  obtained 
from  the  ASW  fusion  cen  ASWOCs  cue  their  MPA  assets  to  investigate  areas 
that  are  likely  to  contain  targets  of  interest.  Remaining  MPAs  also  may  be 
assigned  to  perform  uncued  area  search.  ASWOCs  make  their  MPA  search 
assignments  decisions  regularly  at  user-specified  allocation  intervals.  At  each 
allocation  interval,  the  fusion  center  will  provide  the  ASWOC  with  an  86%  search 
probability  area  (SPA)  for  each  suspected  target  within  the  ASWOCs  ocean  area, 

ASSET  prepares  a  table  that  matches  all  pairs  of  available  MPA  and  SPAs, 
giving  the  times  of  station,  transit  time  to  the  SPA,  and  SPA  size  at  mid-time  on- 
station.  Available  MPAs  include  all  MPAs  on  the  ground  in  a  ready  status  and 
those  MPAs  conducting  uncued  search  that  have  not  yet  reported  a  detection. 

ASSET  first  eliminates  each  SPA/MPA  pairs  whose  projected  SPA  size 
exceeds  a  user-entered  maximum,  and  pairs  whose  computed  time  on-station  is 
less  than  the  user-entered  minimum  (uncued,  divertable  MPA  may  have  a  different 
user-entered  minimum).  From  the  remaining  pairs,  ASSET  selects  the  first 
SPA/MPA  combinations  with  minimum  transit  time.  This  process  continues  until 
either  all  the  SPAs  or  all  the  MPA  have  been  exhausted.  If  any  available  MPA 
remains,  up  to  a  user-entered  maximum  will  be  assigned  to  search  a  user- 
designated  uncued  search  region.  The  MPA's  search  region  is  randomly  located 
inside  this  designated  uncued  area  search  region. 
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[SH  91)  points  out  that  the  way  ASSET  selects  SPA/MPA  pairs  can  be 
improved.  Also  ASSET  only  allocates  MPA  at  the  predetermined  allocation 
intervals. 

C.  IMPROVEMENTS 

[SH  91]  has  suggested  several  alternate  models  of  MPA.  Most  of  the  MPA 
modeling  improvements  mentioned  here  are  derived  from  that  report  but  with 
some  modifications. 

1.  Glimpse  Rate  Model 

All  the  continuous  detection  sensors  in  ASSET  except  the  MPA  use  a 
glimpse  rate  to  determine  detection  opportunities.  An  MPA  Glimpse  Rate  Model 
(MGRM)  approximates  a  continuous-looking  sensor  pattern  that  has  a  probability 
of  detection  (Pd)  of  less  than  1.0  with  a  glimpsing  sensor  region  that  has  a  Pd  of 
1.0  (as  ASSET  currently  does  with  the  tripwire  sensor).  The  sensor  region  would 
then  be  glimpsed  to  provide  a  detection  rate  identical  with  that  obtained  by  a 
continuous  sensor  conducting  a  random  search. 

Detections  for  MGRM  are  based  on  the  random  search  model  where  the 
detection  rate  is  the  ratio  of  the  relative  area  searched  per  unit  time  over  the  total 
area  of  search  region.  Random  search  predicts  that  a  target,  moving  randomly 
through  a  field  of  continuous  stationary  sensors,  will  be  detected  at  a  certain  rate. 
By  glimpsing  the  sensor  field  at  this  detection  rate,  MGRM  can  produce  an 
identical  detection  rate  regardless  of  how  the  target  is  moving. 

a.  Area  of  the  Search  Region 

For  a  cued  search,  the  sonobuoys  will  be  uniformly  placed  through 
the  circular  search  region  of  tracker-correlator's  best  estimate  of  target  position. 
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Otherwise,  the  search  region  will  be  randomly  placed  within  the  user-entered 
uncued  search  region. 

b.  Determination  of  Detection 

When  an  MPA  arrives  on-station,  MGRM  will  calculate  a  glimpse 
rate  and  initial  glimpse  time  for  each  potentially  detectable  target.  As  in  A.SSETs 
other  glimpse  models,  a  target  is  deemed  potentially  detectable  if  its  range  from 
the  MPA  is  less  than  its  maximum  speed  multiplied  by  the  MPA  search  time.  The 
initial  glimpse  time  will  be  determined  from  an  exponential  random  draw  using  the 
inverse  of  the  glimpse  rate  (or  detection  rate)  as  the  mean  glimpse  interval. 
Detection  rate  will  be  computed  as  ASSET  currently  does.  If  the  target  is  within 
the  SPA  when  the  glimpse  occurs  it  is  considered  a  successful  detection. 

c.  Effect  of  False  Contacts 

The  false  contacts  are  generated  by  a  Poisson  process  with  a  user- 
entered  false  contact  rate.  An  exponentially  distributed  random  number  will  be 
drawn  to  determine  the  time  of  the  false  contact.  It  will  transmit  a  false  target 
report  at  the  lime  of  false  alert  if  the  MPA  is  still  on  station  and  has  not  yet  been 
killed.  The  MPA  will  spend  the  same  amount  of  buoys  for  target  classification,  as 
it  does  in  a  normal  engagement.  A  false  contact  will  not  be  processed  if  the  MPA 
is  engaged  at  the  time  of  the  false  alert. 

2.  MPA  Allocation  Model 

a.  Determination  of  SPA/MPA  Matching 

The  allocation  of  the  next  MPA  would  be  to  the  MPA/SPA  pair  with 
the  largest  ratio  of  MPA's  time  on  station  to  SPA  size.  This  will  increase 
cumulative  detection  probability  without  increasing  the  processing  time  to  perform 
an  allocation. 
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b.  A fter  Prosecuting  a  Target 

When  the  MPA  finishes  the  prosecution  of  a  submarine,  it  will  stay 
in  the  search  region  until  the  end  of  its  time  on  station.  This  will  enable  the  MPA 
to  have  opportunities  to  detect  other  targets.  Also  if  time  permits,  the  MPA  can  be 
directed  to  other  cues. 

D,  IMPLEMENTATION  DETAILS 

This  section  provides  implementation  details  for  the  suggested  improvements. 
They  are  categorized  by  their  physical  locations  in  the  application  package 
distributed  by  Metron. 

1.  Structui  e  File 

a.  MPAstruct 

The  MPAstruct  defines  the  operation  characteristics  of  an  MPA. 
One  data  field  is  added: 

engaged?  Indicates  whether  an  MPA  is  currently 

engaging  with  opponent  platform 

2.  Class  ASWOC 

The  ASWOC  class  describes  the  ASWOC's  geographical  location, 
assigned  squadrons,  specific  ocean  area  of  responsibility,  and  operations.  We 
redefine  two  of  its  methods: 

a.  Modified  Functions 

makeMPAassignment :  The  selection  of  next  cue  is  the  cue  with 
the  largest  ratio  of  MPA’s  time  on  station  over  SPA  size. 
getSearchValues  :  SearchValues  now  contains  seaichTime, 
missionCount,  MPA,  SPA  size. 
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3.  Class  MPAsquadron 

The  class  MPAsquadron  defines  variables  of  its  command,  operation 
characteristics,  weapon  status,  and  operations.  We  describe  the  modification  in 
following  sections: 

a.  Modified  Functions 

dispatchMPAtoAraaSaarch  &  diapatchMPAtoCua: 
Schedule  detection  opportunities  for  the  MPA  to  every  possible 
submarine  target  based  on  the  inversion  of  detection  rate  to  the 
target. 

endMPAf  light :  Resets  MPAstruct-engagedp  flag  to  NIL  after 
finishing  its  flight. 

beginSearch :  Schedules  detection  opportunity  for  the  MPA  to 
every  possible  submarine  target  based  on  the  inversion  of  detection 
rate  to  the  target. 

localizeTarget :  Checks  if  the  MPA  is  on  station  and  the 
target  is  active  before  doing  anything. 

loseMPA:  Checks  if  the  MPA  has  not  being  killed  yet  before 
d^Jlng  anything. 
b  New  Functions 

setNextDetectionOpportunity :  Schedules  next  detection 
opportunity  for  MPA. 

processDetectionOpportunity :  Determines  target 
detection;  schedule  next  detection  opportunity. 
targetDetectedp :  Determines  it  the  target  is  detected  by 
checking  if  the  target  is  inside  search  region. 
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targetCovaredp :  Checks  if  the  submarine  is  inside  the  given 
search  region. 

setHextFalseAlarm:  Schedules  false  alarm  event  during 
MPA's  time  on  station  by  making  exponential  draw  from  mean 
time  between  false  alarm. 

procea.'sFalseAlarm:  Transmit  a  false  target  report  if  the 
MPA  is  still  on  station  and  not  being  killed;  expends  buoys  for 
target  classification;  schedule  next  false  alarm  event. 
andEngagament :  Reset  IvIPAstnict-engagedp  flag  to  Nil. 
transmitFalsaAlarmRaport :  Transmit  the  false  target 
report  to  its  command. 

prosecuteTarget ;  Determines  who  wins  during  engagement; 
if  time  permit  and  still  carry  enough  weapons  MPA  will  stay  for 
another  detection  opportunity. 
c.  Deleted  Functions 

getSublnSearchArea 
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V.  CONVERTING  ASSET  TO  CLOS 


flic  ASSET  was  implemented  in  Macintosh  Allegro  Common  Lisp  Version 
1.3.2  which  uses  Object  Lisp  as  object-oriented  extension  of  Common  Lisp. 
Recently  Apple's  support  of  this  version  was  stopped,  and  CLOS  became  the 
standardized  Common  Lisp  package  for  version  2.0  release.  For  the  benefit  of 
future  maintenance  of  the  ASSET,  we  decide  to  convert  the  current  code  to 
CLOS. 

This  chapter  pro''ides  a  summary  of  the  differt  "ices  between  Object  Lisp  and 
CLOS,  the  work  of  converting  current  ASSET  code  to  Cl  .OS,  and  the  results  w'e 
have  achieved. 

A.  DIFFERENCES  BETWEEN  CLOS  AND  OBJECT  LISP 

Machintosh  Common  Lisp  (MCL)  Version  1.3.2  uses  the  object  protocol 
Object  Lisp  which  supports  multiple  mheritance  but  only  simple  method 
combination.  In  Object  Lisp  one  could  mako  an  object  submarine  with  an 
instance  SSN  688,  then  create  SSN  977  which  inherits  from  SSN  688,  and  SSN 
755  which  inherits  form  SSN  977.  CLOS  uses  a  class-instance  protocol.  In 
CLOS  we  cannot  make  a  subclass  from  SSN  688,  since  it  is  not  a  class.  Instead, 
we  must  create  anotiier  class,  possibly  a  subclass  of  submarine. 

Binding  and  scoping  have  changed  substantially.  Instead  of  object  variables, 
there  are  slots  in  CLOS  classes  and  instances.  Rather  than  asking  an  object  to  run 
its  version  of  an  object  function,  CLOS  applies  methods  of  generic  functions  to 
instances.  The  detailed  information  of  difference  between  Object  Lisp  and  CLOS 
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can  be  found  in  the  user's  manual  of  Macintosh  Allegro  Common  Lisp  Version 
1.3.2  and  (ST  90]. 

P.  CONVERTING  TO  CLOS 

1.  Class  Definition 

Generally  in  MCL  1.3,  defob  ject  defines  the  class  hierarchy,  and 
exist  creates  an  object  instance  and  sets  its  own  binding  of  variables  to  initial 
values  by  using  function  have.  Most  frequently,  the  corresponding  material  in 
CLOS  goes  into  the  def class  with  its  slot  specifiers.  The  initialize- 
instance  specifies  the  values  of  slots  that  cannot  initialized  with  initialization 
arguments  or  initial  forms.  It  also  perform  any  other  necessary  initialization.  If 
have  is  used  dynamically  in  the  program,  one  can  transform  uses  of  have  into 
static  slots  and  values  in  defclass. 

2.  Method  Definition 

All  object  functions  must  be  turned  into  generic  functions.  That  is, 
every  piece  of  code  of  the  form 

(defobfun  (Die  submarine)  arguments  body) 
must  be  transformed  to 

(defmethod  Die  (  {submarine  submarine)  r  :r-of-args) 
body) 

3.  Replace  Free  Variables  References 

Since  the  slot- values  and  variables  do  not  access  the  same  namespace,  we 
must  bring  all  free  variable  references  into  the  slot  namespace.  We  have  two 
ways  of  dealing  with  free  variables:  set  them  explicitly  with  slot -value  or  use 
accessor  calls.  Using  accessors  arc  usually  preferable,  since  they  allow  us  to 
change  the  representation  of  our  classes  without  changing  any  of  the  user  code. 
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4.  Remove  Calls  to  Ask 

Instead  of  ask,  we  now  use  slot -value  or  an  access  >t  for  a  value, 
or  a  call  to  a  method.  Depending  on  the  context,  we  should  choose  the 
appropriate  call.  We  must  run  methods  on  instances  of  a  class.  All  calls  must 
now  be  directed  to  instances,  not  to  class  object.  An  Object  Lisp  asking  for  the 
object  function  setForStart 

(ask  submarine  (setForStart)) 

becomes  a  call  to  the  appropriate  method  of  setForStart  for  an  instance  of 
acousticPlatform,  SSN688: 

(setf  SSN688  (make-instance  'submarine)) 

(setForStart  SSN688) 

5.  Create  Instance 

Rather  than  creating  all  new  objects  with  oneof ,  we  create  instances 
with  make-instance.  Similarly,  we  replace  kindof  with  def  class.  Note 
that  in  CLOS,  we  cannot  make  instances  directly  into  classes.  Instead,  we  make 
new  classes  based  on  the  class  we  want  to  specialize,  then  create  instances  of  those 
classes. 

6.  Change  Function  Names 

Numerous  functions  names  and  other  symbols  have  changed  in  both 
Common  Lisp  and  Macintosh  Common  Lisp.  Functions  involving  dialog  items 
or  windows  are  frequently  affected. 

7.  Changes  of  Dialogs 

The  implementation  of  dialog  has  changed  substantially  in  MCL  Version 
2.0.  Dialogs  as  a  separate  class  have  disappeared.  Dialog  items  may  now  be 


added  to  all  views.  Some  functions  have  changed  to  reflect  the  new  definition  of 
dialog. 

C.  RESULTS 

The  source  code  of  ASSEl'  contains  about  29,000  lines.  Most  the  work  of 
converting  it  to  CLOS  can  be  done  by  mechanically  transforming  the  source 
code.  We  used  a  conversion  utility  program  to  help  converting  the  class  and 
method  definitions,  referencing  a  free  variable,  running  a  method,  creating  a 
instance,  and  changing  function  names.  This  utility  program  will  accept  a  text 
file  (of  MACL)  and  try  to  convert  it  to  MCL  program.  Unfortunately,  due  to 
time  constraints  we  are  not  able  to  make  the  utility  program  a  total  solution  for 
the  porting.  To  complete  the  conversion,  we  still  need  to  manually  convert  all 
window  system  to  the  newer  version. 


36 


VI.  COxNCLUSION 


A.  SUMMARY 

In  this  thesis,  we  presented  the  idea  of  object-oriented  simulation.  In 
addition,  we  implemented  the  improvements  suggested  in  previous  evaluations  of 
ASSET  in  NFS.  Some  modifications  to  those  suggestions  were  made  during  the 
implementation  piocess.  Part  of  the  current  ASSET  source  code  is  converted  to 
CLOS  but  more  work  is  needed  to  complete  the  task. 

1.  Detection  Model 

A  compound  Lambda-Sigma  Jump  detection  model  was  implemented  to 
simulate  the  detection  of  submarines.  Also  the  target's  most  detectable  frequency 
and  detection  rate  were  allowed  to  vary  with  environmental  region.  Multiple 
engagements  between  platforms  were  allowed. 

2.  MPA  Model 

A  glimpse  rate  model  was  used  to  determine  detection  opportunities  of 
MPA  and  to  approximate  a  continuous-looking  sensor  pattern.  The  glimpsing 
sensor  field  detects  a  target  which  is  within  the  sensor  region  at  the  time  of  a 
glimpse  with  a  probability  of  detection  of  1.0. 

ASSET  allocates  MPA  to  cues  generated  from  the  tracker-correlator  by 
selecting  SPA/MPA  pairs  using  the  ratio  of  MPA's  time  on-station  to  the  SPA 
size  as  the  selection  criterion.  If  time  permits,  the  MPA  will  stay  on  search 
region  after  prosecuting  a  submarine  for  another  detection  opportunity. 
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3.  Converting  ASSET  code  to  CLOS 

We  investigated  the  feasibility  and  technique  of  converting  current 
ASSET  code  to  CLOS.  We  used  a  conversion  utility  to  help  us  transform  the 
current  ASSET  code  to  CLOS.  Due  to  time  constraints,  only  part  of  the 
converting  task  was  completed. 

B.  FUTURE  WORK 

There  are  several  areas  that  could  be  enhanced  to  improve  the  ASSET.  The 
user  interface  could  be  improved  for  easier  setup  of  a  simulation  architecture.  It 
is  possible  to  speed  up  the  overall  performance  of  the  ASSET  by  porting  it  to  a 
SUN  workstation  with  a  more  powerful  CPU.  Various  modules  such  as  surface 
ASW  platforms,  active  and  nonacoustic  sensors  could  be  added  to  form  a  more 
complete  ASW  simulation.  In  addition,  it  would  be  useful  to  add  machine 
intelligence  features  into  platform  tactics  to  better  represent  smart  platforms. 
The  bottle  neck  of  ASSET  performance  has  not  been  exactly  located.  It  is 
possible  to  further  speed  up  ASSETs  performance. 
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APPENDIX  A  MODIFIED  MODULES  OF  ASSET  CODE 


CHANGE  LOG: 

in  (defstruct  platformAcousticParms) 
new  data  field  added  for  individual  lambda-sigma  process  and 
detection  model: 

(individualSigma  0) 

(lambdaSigmaValue  0) 

(meanlnterval  1000) 

(detectableinterval  1) 

(env-freq-list  nil) 


in  (defstruct  MPAstruct) 
new  data  field  added  for  MPA's  detection  model: 
indicates  whether  a  MPA  is  engaging  with  (or  been  engaged  by)  a 
target 

(cngagedp  nil) 


(defstruct  coast-event 
time 
object 
procedure 
data 

updateList 

) 


sigma  for  individual  fluctuation 

fluctuation  of  target's  signal  excess 
mean  time  between  fluctuation 
target's  detectabl  interval(  1/lambda) 
for  scheduling  detection  event 
environment-frequency  pairs  list  for  a 
target  used  to  determine  its  radiated 
frequjency  in  a  particular  environment 


(defstruct  coast-message 
send-time 
receipt-time 
type 
sender 
content 
size 


;  Time  message  sent. 

;  Time  message  received  at  last  node. 
;  E.g.  'Detect-msg 
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transmission-path 

transmission-count 

) 


(defsmict  obu-contact 
id 

receipt-time 

track-association 

sensor 

categorization 

spadal'data 

altitude 

HFDFp 

) 


(defstruct  obu-track 
id 

number-of-contacts 

head-state-covariance 

head-contact-id 

head-spatial-data 

taii-state-covariance 

tail-contact-id 

tail-spatial-data 

altitude 

(contacts  nil  :type  list) 

) 


(defstruct  obu-data-field 
type 

obs-time 

lat 

Ing 

brg 

maj 

min 

cse 
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cse-unc 

spJ 

spd-unc 

) 


(defstruct  obu-state-field 
ref-lat 
ref-lng 

(state  nil  ;type  (array  float  4  )) 
(covariance  nil  :type  (array  float  '(4  4))) 

) 


(defstruct  obu-sensor 
name 

single-report-to-trackp 

) 


(defstruct  obu-category 
pinnedp 
lockedp 
assoc-count 
cluster-count 
ambiguity-list 
ambiguity-resolution 
ambiguity-resolution-time 

) 


(defstruct  sen-detect-parms 
sensorName 
ga 

bw 

pd 

pfa 

rho 

msize 
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min-delay 

max-delay 

type 

pos-unc 

spd-unc 

cse-unc 

) 


(defstnict  target-state-parms 
targetName 
lat 
Ing 
cse 
spd 

detec  tablep 
communicatingp 

) 


(defstnict  platformStateParms 
(platformName  "") 

(platformType  ’general) 

(classType  nil) 

(gcncr^CommTypc  ’surface) 
(targetType  nil) 

(side  nil) 

(opponents  nil) 

(lat  0) 

(Ing  0) 

(hgt  0) 

(cse  0) 

(spd  0) 

(nav  0) 

(activep  t) 

(communicatingp  nil) 

(engagedp  nil) 

(engagementinterval  nil) 

(engagedPlatform  nil) 

(inTrailp  nil) 

(trailedPlatform  nil) 
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(cuedp  nil) 

(cueinginterval  nil) 
(interceptingp  nil) 

(targetObject  nil) 

(postEngageMotionOrder  nil) 
(deletedp  nil) 


) 


(defstruct  platformAcousticParms 
(sigFreq  nil) 

(spdThresh  0) 

(RNslow  0) 

(RNfast  0) 

(SNlist  nil) 

(DI-RD-list  nil) 

(sigFreq  nil) 


) 


new  data  field  added  for  individual  lambda-sigma 
process 


(individualSigma  0) 
(lambdaSigmaValue  0) 
(meaninterval  1000) 
(detectableinterval  1) 

(env-freq-list  nil) 


ifluctuation  of  target’s  signal  excess 
;mean  time  between  fluctuation 
itarget's  detectabl  interval(l /lambda) 
;for  scheduling  detection  event 
;environment-frequency  pairs  list  for  a 
;target;  used  to  determine  its  radiated 
;frequjency  in  a  particular  environment 


(defstruct  locationOrder 
lat 
Ing 
hgt 

) 
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(defstruct  movementOrder 
time 
cse 
tLat 
tLng 
spd 
nav 
area 

) 


(defstruct  patrolOrder 
time 
region 

patrolDuration 

) 


(defstruct  transiiOrder 
time 
track 

legnumber 

) 


(defstruct  planOrder 
time 
plan 

stepnumber 

duration 

legnumber 

repositionp 

) 


(defstruct  region-parms 
regionName 
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cPoints 

minMax 

minLat 

maxLat 

minLng 

maxLng 

cenLat 

cenLng 

gridpts 

ptArea 

total  Area 

checkpts 

thePoies 

theNormals 

theRngs 

theBrgs 

cVectors 

sVectors  ) 


(defstruct  posit 
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lat 

Ing 

hgt 

time 


(defstnict  MPAstruct 
ASWOC 
squadron 
baseLat 
baseLng 

(readyp  t) 

(searchingp  nil) 

searchStartTime 

searchEndTime 

missionEndTime 

lat 

Ing 

radius 

searchRegion 
(buoyCount  0) 
(torpedoCount  0) 

(totalFlightHours  0) 
(failureTime  0) 

(missionCount  0) 


new  data  field  added  for  MPA's  detection  model; 

;;  indicates  whether  a  MPA  is  engaging  with(or  been  engaged 
;;  by)  a  target 


(engagedp  nil) 

drawnMission 

keepOutOrder 

nextContact 

kiiledp 


(defstruct  spaceSensorStruct 
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keyl 

key2 

biioleanKeys 

keyl  bins 

key2bins 

velocityBiiis 

deplhBins 

probabilities 

) 


(defstruct  se'isorPd 
FAR 
WirePd 
buoyPdVec 
boaiPd  Array 

) 


(defstruct  NAArea 
NA-index 
windSpeed 
windDirection 
visibility 
surfAirTemp 
seaSurfTemp 
cloudCover 
JerlovType 
bioluminesence 

) 


(defstruct  forceStruct 
side 

sideColor 

menuTitle 

forceDefinition 

workSpaceList 

workSpaceFile 

classFile 
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noiseFrequencyFiL 

noi.seFrequencir>’ 

objectKeyList 

instanceList 

opponentList 

linkLis'. 

connectList 

) 


(defStruct  objectStruct 
template 
index 

(classNameList  ()) 

(classList  0) 

(objectNameList  ()) 

(objectList  ()) 

menuTiile 

disabledp 

dialogText 

classObject 

platformType 

initList 

) 


(defstruct  subobject-category 
name 

parent-subobject 

descendent-subobject-list 

) 


(defstruct  radio-button 
symbol 
value 

) 
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;;  (4)  COAST  MOE  CALCULATION  STRUCTURES: 


(defstruct  coast-moe 

(number-of-data-elements  0) 

(sum-of'data  0) 

(sum-of'Squared-data  0) 

(data  nil) 

(abscissa  0) 

) 


(defstruct  MOEstruct 

menuTitle 

headings 

CGheadings 

side 

eventType 

initialcount 

count 

MOEarray 

cuirentlndex 

elapsedTime 

) 


(defstruct  MOEeventStruct 
repNumber 
time 
eventID 
object  1 
object2 
data 

) 


(defstruct  MOEobjectStruct 
side 
type 
class 
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name 


CHANGE  LOG: 


lambdaSigmaProcess-rb  added 
radio  button  for  inputting  lamda-sigma  parameters 

env-freq-rb  added 

radion  button  for  inputting  environment-frequency  pairs  of 
a  aconsticplatform 


;; . LOAD  "Object  Editor  RBs"  Folder  Contents; 

(let  ((directory  "coastLibrary;FolderH  -  RBs:") 

(files  (list  "general  RB  type" 

"data-edit-rb" 

"Lat  Lng  RB" 

"Select  and  Modify  RB" 

"Table  Selector  RB" 

"Build  Connectivities  RB")) 

( dolist  (i  files) 

(load  (make-pathname  :directory  directory  :name  i)) 

) 


(ask  objectEditor 

(have  'rb-installation-list  nil) 


(have  'Sensor-rb 

(oneof  select-&-modify-rb 

:dialog-item-text  "Sensors" 

:window-title  "Select  Sensors" 
;dialog-item-position  (make-point  -1  391) 
:dialog-item-size  (make-point  140  16) 
.radio-button-pushed-p  nil 
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:radio-button-c  luster  '0 
:rb-  symbol  'Sensor-rb) 

) 

(have  'Link-rb 

(oneof  select-&-modify-rb 

:dialog-item-icxt  "C3  Links” 
:window-titIe  "Select  Links" 
:dialog-item-position  (make-point  -1  408) 
:dialog-item-size  (make-point  140  16) 
;radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:rb-symbol  'Link-rb) 

) 


(have  'Connect-rb 

(oneof  ConnectivityRadioButton 

:dialog-item-text  "C3  Connectivities" 
:dialog*item-position  (make-point  -1  426) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:rb-symbol  'Connect-rb) 

) 


(have  'Location-rb 
(oneof  lat-lon-rb 

:dialog-item-text  "Location" 
:dialog-item'position  (make-point  -1  338) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
;rb-symbol  'Location-rb) 

) 


(have  'Region-rb 

(oneof  table-selector-rb 
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:dialog-item-text  "Region" 

:dialog-item-position  (make-point  -1  374) 
;dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nit 
:radio-button-cluster  '0 
:rb-symbol  'Region-rb 

:table-data-form  '(ask  regionManager  rgnList)) 

) 

(have  'Tripwire-rb 
(oneof  table-selector-rb 

:dialog-item-text  "Tripwire" 
rdialog-item-posidon  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
;radio-button-pushed-p  nil 
:radio-button-cluster  ‘0 
:rb-symbol  'Tripwire-rb 
:tabIe-data-form  '(ask  trackManager  trackList)) 

) 

(have  'Sensor-Region-rb 
(oneof  lable-selector-rb 

;dialog-item-texi  "Sensor  Region" 
:dialog-item-position  (make-point  -1  374) 
;dia]og-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:rb-symbol  'Sensor-Region-rb 
:table-data-form  (ask  regionManager  rgnList)) 

) 

(have  'Operaiing-Region-rb 
(oneof  table-selector-rb 

:dialog-item-text  "Operating  Region" 
;dialog-item-position  (make-point  -1  374) 
;dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  0 
:rb-symbol  'Operating-Region-rb 
rtable-data-form  '(ask  regionManager  rgnList)) 

) 

(have  'Area-Search-Region-rb 
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(oneof  table-selector-rb 

:dialog-itenrt-text  "Area  Search  Region" 
:dialog-iiem-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:rb-symbol  'Area-Search-Region-rb 
:table-data-foTm  '(ask  regionManager  rgnList)) 

) 

(have  'Plan-rb 

(oneof  table-selector-rb 

:dialog-itcm-tcxt  "Motion  Plan" 

:dialog-item-position  (make-point  -1  374) 
rdialog-item-size  (make-|>oint  140  16) 
:radio-button-pushed-p  nil 
:radio-buiion-cluster  '0 
:rb-symbol  'Plan-rb 

:table-data-form  '(ask  PlanManager  motionPlanList)) 

) 

(have  'Replenishment-rb 
(oneof  table-selector-rb 

:dialog-item-texi  "Replenishment  Plan" 

;dialog-item-position  (make-point  -1  374) 

rdialog-item-size  (make-point  140  16) 

:radio-button-pushed-p  nil 

rradio-button-cluster  '0 

:rb-symbol  'Replenishment-rb 

rtable-dala-form  '(ask  PlanManager  motionPlanList)) 

) 

(have  'commander-rb 
(oneof  table-selector-rb 

rdialog-item-text  "Commander" 
rdialog-item-position  (make-point  -1  374) 
rdialog-item-size  (make-point  140  16) 
rradio-button-pushed-p  nil 
rradio-button-cluster  '0 
rrb-symbol  'commander-rb 
rtable-data-form 
'(let  ((side  (ask  my-dialog 

(ask  currentObject 
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) 


(platformStatePanns-side  platfoimState)))) 
(indexList  (list  (ask  my-dialog  (ask  currentObject 

commandlndex))))) 

(getNamesFromlndices  side  indexList))) 


(have  'fusionCenter-rb 
(oneof  table-selector-rb 

:diaIog-item-text  "Fusion  Center" 

:dialog>item-position  (make-point -1  374) 

rdialog-item-size  (make-point  140  16) 

:radio-button-pushed-p  nil 

:radio-button-cluster  ‘0 

:rb-symbol  'fusionCenter-rb 

:tabIe-data-form 

'(let  ((side  (ask  my-dialog 

(ask  currentObject 

(platformStateParms-side  plalformState)))) 
(indexList  (list  (ask  my-dialog  (ask  currentObject 

fusionindex))))) 

(getNamesFromlndices  side  indexList))) 

) 


(have  'frequency-rb 
(oneof  table-selector-rb 

:dialog-iiem-text  "Signature  Freq" 
:dialog-item-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-butfon-cluster  '0 
:rb-symbol  'frequency-rb 
:table-daia-form 
'(ask  my-dialog 
(ask  currentObject 

(getFreqList  (platformStateParms-side 

platformState))))) 


) 


(have  'subMissions-rb 
(oneof  table-selector-rb 

:dialog-item-text  "Select  Mission" 
:dialog-item-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
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) 


:radio-button-pushed-p  nil 

rradio-button'Cluster  '0 

;rb-syinbol  'subMissions-rb 

:table*data-form  '(^k  my-dialog  (ask  currentObject 

subMissionList))) 


(have  'AcousSig*rb 
(oneof  data-edit-rb 

:dialog-item*text  "Acoustic  Sig" 
:dialog-item-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nit 
rradio-button-cluster  '0 
:data-key  :acousticSig 

) 

) 


♦  ♦ 

;;  new  radio  button  for  acoustic  platform 


(have  'LambdaSigmaProcess-rb 
(oneof  data-edit-rb 

:dialog-itcm-tcxt  "Lambda-Sigma  Process" 
:dialog-item-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
;data-key  :LambdaSigmaProcess 

) 

) 


(have  'SubmarineWeaponsl-rb 
(oneof  data-edit-rb 

:dia]og-item-text  "Submarine  Weapons  1" 
:dialog-item-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 


56 


:radio-button-cluster  ’0 
:data-kty  :SubmarineWeaponsl 

) 

) 

(have  'SubmarineWeapons2-rb 
(oneof  data-edit-rb 

:dialog-item*text  "Submarine  Weapons  2" 
:dialog-item-position  (make-point  -1  374) 
;dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-butlon-cluster  '0 
:data-key  :SubmarineWeapons2 

) 

) 


(have  'SubmarineTactics-rb 
(oneof  data-edit-rb 

:dialog-item-text  "Submarine  Tactics" 
:dialog-item-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-buiton-pushed-p  nil 
:radio-button-cluster  '0 
:data-key  .SubmarineTactics 


) 


) 


(have  'sensorParms-rb 
(oneof  data-edit-rb 

:dialog-item-text  "Detection  Parameters" 
rdialog-item-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:data-key  :sensorParms 


) 


) 


(have  'MPASquadron-rb 
(oneof  data-^it-rb 

:dialog-item-text  "Squadron  Parameters" 
rdialog-item-position  (make-point  -1  374) 
rdialog-item-sizc  (make-point  140  16) 
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:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:data-key  :MPASquadron 


) 

) 

(have  'MPAstores-rb 
(oneof  data-edit-rb 

:dialog-iiem-text  "Stores  Parameters" 
:dialog-item-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:data-key  :MPAstores 

) 

) 


(have  'general-motion-rb 
(oneof  data-edit-rb 

:dialog-item-text  "General  Motion" 
:dialog-item-position  (make-point  -1  374) 
:dialog-item-sizc  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:data-key  :generalMotion 

) 

) 

(have  'reports-rb 
(oneof  data-edit-rb 

:dialog-item-text  "Detection  Reports" 
:dialog-item-position  (make-point  -1  374) 
;dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:data-key  :reports 

) 

) 

(have  'orbitParms-rb 
(oneof  data-edit-rb 

;dialog-item-text  "Orbital  Parameters" 
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) 

) 


rdialog-itenri-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:data-key  rorbitParms 


(have  'SOAallocation-rb 
(oneof  data-edit-rb 

idialog-item-text  "Alloc  Parameters" 
:dialog-item-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
;radio-button-cluster  *0 
;data-key  rSOAallocationParms 


) 


) 


(have  'ASWOCallocation-rb 
(oneof  data-edit-rb 

:dialog-item-text  "Alloc  Parameters" 
:dialog-itcm-position  (make-point  -1  374) 
:diaIog-item-size  (make-point  140  16) 
:radio-button-pu3hed-p  nil 
:radio-button-cIuster  '0 
:data-key  ;ASWOCallocationParms 


) 


) 


(have  'HFDF-rb 
(oneof  data-edit-rb 

:dialog-item-text  "HFDF  Detection" 
:dialog-item-position  (make-point  -1  374) 
:dia]og-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
;radio-button-cluster  '0 
:data-key  :HFDFparms 


) 


) 


(have  'minefield-rb 
(oneof  data-edit-rb 
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:diaIog-item-text  "Minefield  ParanKters" 
:dialog-item-position  (make-point  -1  374) 
'.dialog-item- size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:data-key  :MineFieldParms 


) 

(have  'tripwirePanns-rb 
(oneof  data-edit-rb 

:dialog-item-text  '  Tripwire  Parameters" 
:dialog-item-position  (make-point  - 1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  ‘0 
:data-key  :TripwirePanns 

) 

) 


(have  'submarineSensorPD-rb 
(oneof  coast-RB-type 

:dialog-item-text  "Sub  Sensor  PDs" 

:dialog-item-position  (make-point  - 1  374) 
;dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:dialog-item-action 

’(ask  (ask  my-dialog  currentObject)  (editSubmarinePDs)) 

) 

) 

(have  'surfaceSensorPD-rb 
(oneof  coast-RB-type 

:dialog-item-iext  "Surf  Sensor  PDs" 

:dialog-itern-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:dialog-item-action 

'(ask  (ask  my-dialog  currentObject)  (editSurfacePDs)) 
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(have  'subSubPk*rb 
(oneof  coasl-RB-type 

:dialog-item-text  "Sub  Sub  PKs" 

:dialog-item-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radiO'button-pushed-p  nil 
:radio-buiton-cluster  '0 
;dialog-item-action 

'(ask  (ask  my-dialog  currentObject)  (editSubSubPks)) 


(have  'MPASubPk-rb 
(oneof  coast-RB-type 

:dialog-item-text  "MPA  Sub  PK-" 

:diaIog-item-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
:radio-button-pushed  t  nil 
:radio-bulton-cluster  C 
:dialog-item'action 

'(ask  (ask  my-dialog  currentObject)  (editMPASutPks)) 


(have  minefie!dRange-rb 
(oneof  coast-RB-type 

:dialog-item-text  "Sub  Mine  Ranges" 

;dialog-item-position  (make-point  -1  374) 

;dialog-item-size  (make-point  140  16) 
.radio-button-pushed-p  nil 
:radio-buiton-cluster  0 
:dialog-item-action 

'(ask  (ask  my-dialog  currentObject)  (editMincFieldRanges)) 


(have  'minefieldPk-rb 
(oneof  coast-RB-type 

:dialog-item-text  "Minefield  Pks" 
.dialog-item-position  (make-point  -1  374) 
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) 


) 


:dialog-item-size  (make-point  140  16) 
:radio-faLtton-pushed-p  nil 
:radio-button-cluster  '0 
:dialog-item-action 

'(ask  (ask  my-dialog  currentObject)  (editMinePks)) 


(have  'SN-rb 
(oneof  coast-RB-type 

;dialog-item-text  "Self  Noise" 
:dialog-item-position  (make-point  -1  374) 
:dialog-item-size  (make-point  140  16) 
;radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:dialog-item-action 

'(ask  (ask  my-dialog  currentObject)  (editSN)) 

) 


) 

(have  'DI-RD-rb 
(or.eof  coast-RB-type 

:dialog-item-text  "DI,  RD" 

;dialog-item-position  (make-point  -1  374) 
:diaiog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
:radio-button-cluster  '0 
:dialog-item-action 

'(ask  (ask  my-dialog  currentObject)  (editDI-RD)) 


) 


) 


;;  new  radio  button  for  acoustic  platform 


(have  'ENV-FREQ-rb 
(oneof  coast-RB-type 

:dialog-item-text  "Freq.  -  Environment" 
:dialog-item-position  (make-point  -1  374) 
;diaJog-item-size  (make-point  140  16) 
:radio-button-pushed-p  nil 
.'radio-button-cluster  '0 
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:dialog-iiem-action 

'(ask  (ask  my-dialog  currentObject)  (editEnv-Freq)) 


) 

) 
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CHANGE LOG: 


data  field  "sigma"  in  EXIST  function  deleted; 
getPL-AN-Sigma  modified 

target's  radiated  frequency  is  used  on  deciding  PL  &  AN 
getSweepWidth  modified 

target's  radiated  frequency  is  used  on  deciding  PL  &  AN 
which  in  turn  decide  sweep  width 


(setq  environmentManager  (kindof  nil)) 


(defobfun  (exist  environmentManager)(init-Iist) 

(usual-exist  init-list) 

(have  'allPLfrequencies  nil) 

(have  'EnvRegionList  nil)  ;form  is:  ({  Env-rgn  objl  )  {  Env-rgn 

;obj2  }  ...) 

(have  'PL-list  nil) 

;PL-list  will  be  in  following  form: 

;  (Freql  (PL-11  PL-12  .,.)  Freq2  (PL-21  PL- 22  ...)  .  .  .) 

;  where  PL-ij  is  in  form:  C'PLname"  "PLprimaryRgn"  PL-array 
;  FreqAtWhichPLdefined) 

;  NOTE:  FreqAtWhichPLdefined  will  equal  FreqK  in  PL-lisi 
immediately  proceeding 

;  the  proploss  list  (ie  all  PL-ij's  in  (getf  PL-list  FreqK)  will 

;  be:  ("PLname"  "PLprimaryRgn"  PL-array  FreqK)  ). 

:  This  info  s  redundant,  but  necessary  since  the  proploss, 

PL-ij,  can  be  assigned  to  an  env  rgn  at  *any*  frequency, 

:  Thus,  in  order  to  re-constitute  the  proploss  we  need  to 

:  know  the  frequency  at  which  it  was  defined  (not  the 

;  frequency  at  which  it  is  being  used). 

(have  'defaultAN  60,0) 
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(have  'defaultlnitialSphericalPL  66.2) 

(have  'defaultSpherical  20.0) 

(have  'defaullInitialCylindricaIPL  31.1) 

(have  'defaultCylindrical  10.0) 

(have  'envFile  (make-pathname 


:directory  "coastData;Environment  Folder:" 
:name  "EnvironmentFile")) 

(printStatus  "LOADING  ENVIRONMENTAL  DATA"  2) 
(readEnvManager) 


) 


(defobfun  (writeEnvManager  environmentManager)  () 
(wiih-open-file  (toStream  envFile 
;direction  loutput 
:if-exists  rsupersede 
lif-does-not-exist  :create 

) 

(fresh-line  toStream) 

(write  allPLfrequencies  :stream  toStream) 

(terpri  toStream) 

(write  defaultAN  ;stream  toStream) 

(terpri  toStream) 

(write  defaultInitialSphericaPL  :stream  toStream) 
(terpri  toStream) 

(write  defaultSpherical  :stream  toStre<un) 

(terpri  toStream) 

(write  defaultCylindrical  :.stream  toStream) 

(terpri  toStream) 

(write  sigma  :stream  toStream) 

(terpri  toStream) 

(write  PL-list  :stream  toStream  :escape  t 
:level  nil  :length  nil  rarray  t) 

(terpri  toStream) 

;Save  #  of  Env  Rgns  in  file: 

(write  (list-length  EnvRegionList)  :stream  toStream) 
(dolist  (thisEnvRgn  EnvRegionList) 

(ask  thisEnvRgn  (writeEnvRgn  toStream)) 

) 

(terpri  toStream))) 
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(defobfun  (readEnvManager  environmentManager)  () 

(unless  (probe-file  envFile) 

(print  "Warning  from  (readEnvManager  environmentManager);  file 

does  not  exist;") 

(print  envFile) 

(rctum-from  readEnvManager  nil) 

) 

(with-open-file  (fromStream  envFile 
;direction  :input 
:if-does-not-cxist  nil 

) 

(let  ((flength  (file-length  fromStream)) 

(envManager  (selO) 
numberOfRgns  thisRgn 
) 

(if  (or  (null  flength)  (<  flength  5)) 

(returri-from  readEnvManager  nil) 

) 

(setq  allPLfrequencies  (read  fromStream  nil  nil)) 

(setq  defaultAN  (read  fromStream  nil  nil)) 

(setq  defaultInitialSphericaPL  (read  fromStream  nil  nil)) 

(setq  defaultSpherical  (read  fromStream  nil  nil)) 

(setq  defaultCylindrical  (read  fromStream  nil  nil)) 

(setq  sigma  (read  fromStream  nil  nil)) 

(setq  PL-list  (read  fromStream  nil  nil)) 

(setq  numberOfRgns  (read  fromStream  nil  nil)) 

(setq  EnvRegionList  nil) 

(dotimes  (k  numberOfRgns) 

(setq  thisRgn  (oneof  EnvRegionType)) 

(ask  thisRgn  (readEnvRgn  fromStream  envManager)) 

(setq  EnvRegionList  (cons  thisRgn  EnvRegionList)) 

) 

(setq  EnvRegionList  (reverse  EnvRegionList)) 

) 

) 

) 


modified  function 
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(defobfun  (getPL-AN-Sigma  environmentManager)(searcherObj  targetObj) 
(let*  ((sPos  (ask  scarcherObj 

(list  (platformStateParms-lat  platformState) 
(platformStatcParms-Ing  platformState)))) 

(tPos  (ask  targetObj 

(list  (platforms  tateParms-lat  platformState) 
(platformStateParms-lng  platformState)))) 

(mg  (car  (getRngBrg  (car  sPos)(cadr  sPos) 

(car  tPos)(cadr  tPos)  *greatCircle*))) 


new  local  data  used  to  obtain  target's  most 
detectable  radiated  frequency 


(targetEnvRgn 

(tagetEnvRgnName 

(searcherEnvRgn 

(ambNoise 

(PL-val 

(AN-PL-pair 

) 


nil) 

nil) 

nil) 

nil) 

nil) 

nil)freq  ;CHANGE 


;;  this  block  modified  to  incorporat  target's  radiated  frequency 
on  deciding  PL  &  AN 


(dolist  (thisEnvRgn  EnvRegionList) 

(when  (objectInRegionp  scarcherObj  (ask  thisEnvRgn  geoRegion)) 
(setq  searcherEnvRgn  thisEnvRgn) 

(return) 

) 

) 

(when  searcherEnvRgn 
(dolist  (thisEnvRgn  EnvRegionList) 

(when  (objectInRegionp  targetobj  (ask  thisEnvRgn  geoRegion)) 
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(setq  targeiEnvRgn  ihisEnvRgn) 

(return) 

)) 

(setq  targetEnvRgnName  (ask  targetEnvRgn  envRgnName)) 

(setq  freq  (first  (gctf  (ask  targctobj 
(platformAcousticParms-env-freq-list 

platformAcoustics))  targetEnvRgnName))) 

(setq  AN-PL-pair  (ask  searcherEnvRgn  (getf  ambNoise&l^opLoss 

freq))) 

(setq  ambNoise  (first  AN-PL-pair)) 

(if  (setq  PL-array  (third  (second  AN-PL-pair))) 

(setq  PL-val  (aref  PL-array  (round  mg))) 

) 

) 


(if  (null  ambNoise)  (setq  ambNoise  defaultAN)) 

(if  (null  PL-val)  (setq  PL-val  (getDefaultSpherical  mg))) 
(retum-from  getPL-AN-Sigma 

(values  PL-val  ambNoise  sigma)) 

) 

) 


modified  function 


(defobfun  (getSweepWidth  environmentManager)(target  NP) 
(let*  ( 


;;  new  local  data  used  to  obtain  target's  most 
;;  detectable  radiated  frequency 


(targetEnvRgn  nil) 
(tagetEnvRgnName  nil) 
(freq  nil) 

(AN-PL-pair  nil) 
(ambNoise  nil) 
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(PL-array  nil) 
(FOM  nil) 
(detectionRange  0) 
) 


»  * 

;;  this  block  modified  to  incorporat  target's 
;;  radiated  frequency  on  deciding  sweep  width 


(dolist  (thisEnvRgn  EnvRegionList) 

(when  (objectInRegionp  target  (ask  diisEnvRgn  geoRegion)) 

(setq  targetEnvRgn  thisEnvRgn) 

(return) 

) 

) 

(when  targetEnvRgn 

(setq  targetEnvRgnName  (ask  targetEnvRgn  envRgnName)) 

(setq  freq  (first  (getf  (ask  target 

(platformAcousticParms-env-freq-list 

platformAcoustics)) 

targetEnvRgnName))) 

(setq  AN-PL-pair  (ask  targetEnvRgn 

(getf  ambNoise&PropLoss  freq))) 

(setq  ambNoise  (first  AN-PL-pair)) 

(setq  PL-array  (third  (second  AN-PL-pair))) 


(if  (null  ambNoise)  (setq  ambNoise  defauliAN)) 

(setq  FOM  (-  NP  ambNoise)) 

(if  (null  PL-array) 

(setq  detectionRange  (RngFromPL&Spreading  FOM  2)) 
(dotimes  (i  (car  (array-dimensions  PL-array))) 

(if  (aref  PL-array  i) 

(if  (<  (aref  PL-array  i)  FOM) 

(setq  detectionRange  i)) 

) 

) 

) 

(retum-from  getSweepWidth 
(*  2  detectionRange)) 
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) 


(defobfun  (getDefauItSpherical  environmentManagerXmg) 

(if  (<  mg  0.001  )(retum-from  getDefauItSpherical  defaultlnitialSphericalPL)) 
(retum-from  getDefauItSpherical 

(+  defaultlnitialSphericalPL  (*  defaultSpherical  (log  mg  10))) 

) 

) 


(defobfun  (RngFromPL&Spreading  environmentManager)  (PropLoss 
SpreadExpt) 

;;  The  spreading  law  exponent  is  1  for  cylindrical,  2  for  Spherical 
The  PL=0  distance  is  assumed  to  be  one  yard,  Inm  =  2027  yds. 

;;  Returns  range  in  nm 
(if  (<  PropLoss  0) 

(retum-from  RngFromPL&Spreading  0.001) 

(retum-from  RngFromPL&Spreading 

(float  (/  (expt  10  (/  PropLoss  (*  10  SpreadExpt))) 

2027)) 

) 

) 

) 


(defobfun  (geiPLnames  environmentManager)  (freq) 

(let  ((PLsAtFreq  (getf  PL-list  freq  "Not  Found")) 

) 

(when  (equalp  PLsAtFreq  "Not  Found") 

(princ  "Warning  from  getPLnames:  no  PLs  at  given  freq  ") 
(princ  freq)  (terpri) 

(retum-from  getPLnames  nil) 

) 

(retum-from  getPLnames  (mapcar  #Tirst  PLsAtFreq)) 
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(defobfun  (getPLregions  environmentManager)  (freq) 

(let  ((PLsAtFreq  (getf  PL-list  freq  "Not  Found")) 

) 

(when  (equalp  PLsAtFreq  "Not  Found") 

(princ  "Warning  from  getPLregions:  no  PLs  at  given  freq  ") 
(princ  freq)  (terpri) 

(retum-from  getPLregions  nil) 

) 

(retum-from  getPLregions  (mapcar  #'second  PLsAtFreq)) 

) 

) 


(defobfun  (getPLarray  environmentManager)  (freq  PLname) 

(let  ((PLsAtFreq  (getf  PL-list  freq)) 

) 

(dolist  (PL  PLsAtFreq) 

(if  (string-equal  (first  PL)  PLname) 

(retum-from  getPLarray  (third  PL))) 

) 

(princ  "Warning  from  getPLarray:  no  PL  with  freq,  name:  ") 
(print  (list  freq  PLname))  (terpri) 

(retum-from  getPLarray  nil) 

) 

) 


(defobfun  (geiPL  environmentManager)  (freq  PLname) 

(let  ((PLsAtFreq  (getf  PL-list  freq)) 

) 

(dolist  (PL  PLsAtFreq) 

(if  (string-equal  (first  PL)  PLname) 

(retum-from  getPL  PL)) 

) 

(princ  "Warning  from  getPL:  no  PL  with  freq,  name:  ") 
(print  (list  freq  PLname))  (terpri) 


(retum-from  geiPL  nil) 

) 

) 


(defobfun  (putPL  environmeniManager)  (PL  freq) 
(let  ((PLsAtFreq  (getf  PL-list  freq  "Not  Found")) 

) 

(when  (equalp  PLsAtFreq  "Not  Found") 

(princ  "Warning  from  putPL;  bad  freq  =  ") 
(princ  freq)  (terpri) 

(retum-from  putPL  nil) 

) 

(setf  (getf  PL-list  freq)  (cons  PL  PLsAtFreq)) 

) 

) 
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;•  CHANGE  LOG: 

»> 

;;  new  data  fields  and  input  dialogs  of  globallambda-sigma  jump  model 
;;  for  signal  excess  calculation  added;  user  is  required  to  specify 
;;  them  during  data  input  stage: 

;;  globalLambdaSigma Value 
;;  globalSigma 
;;  globalMeanlnterval 
1 1 

;;  setUpForReplication  modified 

initializes  globalLambdaSigmaValue  and  schedules  event  in 
exponentially  separate  time  to  update  its  value 

*♦ 

;;  updateGlobalLambdaSigmaValue  added  to  update  global  lambda 
;;  sigma  value  in  exponentially  separate  time 

«  9 

;;  FOMglimpseInterval  and  its  input  dialog  deleted; 

;;  now  using  target's  Deiectableinterval  to  schedule  next  detection 
;;  event 

99 

;;  setNextDetectionOpportunity  modified 

;;  it  now  allows  engaged  target  being  detectable  to  searcher  that  is 
;;  not  engaged:  also  check  range  when  scheduling  next  detection 
;;  opportunity  for  engaged  platforms 

99 

;;  processDetectionOpportunity  modified 

;;  engaged  searchers  searcher  won't  detect  other  targets; 

;;  engaged  targets  is  detectable  to  other  searchers 

99 

;;  getDetectionOppResult  modified 

;;  determines  target's  most  detectable  radiated  frequecy; 

;;  using  global  and  individual  lambda-sigma-value  to  adjust 
;;  MSEfmean  signal  excess) 

99 

;;  processPlatformKill  modified 
;:  check  if  it  was  killed  before 
♦ » 

;;  setNextFalseAlarm  and  processFalse Alarm  modified 
;;  check  if  the  searcher  was  killed  before  doing  anything 
;;  it  won’t  set  the  engaged  flag  of  searcher;  so  it  can  deal 
;;  with  real  target  while  it  is  processing  FA 
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(setq  Umpire  (oneof  nil)) 


(defobfun  (exist  Umpire)  (init-list) 
(usual-exist  init-list) 

(have  prehostilitiesDuration  50) 

(have  hostilitiesDuration  200) 


(have  'globalSigma  0) 

(have  'globalMeanInterval  1000)  ;mean  time  between  change 

;of  lambda-sigma  signal  excess 
(have  'globalLambdaSigmaValue  0)  ;lambda-sigma  signal 

•.excess 


f***»»*mt»*i»»»i**»i*»»»*i*»*»*»i*»»*t»t»t*i***»**»*tt»**»»>»»**»»****»»*»»*»*»* 


(have  'maxClosingSpeed  50) 

(have  'maxSubDetectionRange  100) 

(have  'maxSurfaceDetectionRange  100) 

(have  'hostilitiesp  nil) 


(have  ’name  "Simulation  Umpire  ") 

(have  'btitle  "EDIT  Umpire  Parameters") 

(have  etitle  "EDIT  Umpire  Parameters") 

(have  daiaNameList 

'(:main  (prehostilitiesDuration  hostilitiesDuration 
globalSigma  globalMeanInterval))) 

(have  dataValueList 

'(:main  (.prehostilitiesDuration  .hostilitiesDuration 
.globalSigma  .globalMeanInterval))) 

(have  'dataTypeList 

'(:main  (data  data  data  data))) 

(have  'dataTextList 

'(:main  ("Prehostilities  Duration"  "Hostilities  Duration" 

"Global  Sigma"  "Global  Mean  Int"))) 

(have  dataTemplateText 

'(:inain  (("Hours"  "10  -  unlimited] ")("Hours"  "[0  -  unlinuted]") 
("decibles"  "[0  -  20J")  ("Hours  (1/Lambda)" 

"(0.25  -  unlimited]")))) 
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(have  'dummyRBdata  nil) 

(have  ’Rb_daianamelist  nil) 

(have  ’Rb_datavaluelist  nil) 

(have  'Rb.designatorlist  nil) 

(have  'UmpireDatafile 
(make-pathname 

idireciory  "coastData;Umpire  Folder:" 
:name  "Umpire  Data  File")) 
(printSiatus  "LOADING  UMPIRE  DATA"  2) 
(recallUmpireData) 

) 


(defobfun  (editParameters  Umpire)  () 
(let  ((mySelf  (selO)) 

(editObjeci  mySelf  :main) 

) 

) 


(defobfun  (recallUmpireData  Umpire)  () 
(with'Open-file  (rstream  umpireDataFile 
:direciion  .input 
:if-does-not-exist  :create) 

(setq  dataValueList  (read  rstream  nil  dataValueList)) 
(mapc  #'(lambda(x  y) 

(if  (listp  x)(mapc  #'set  x  y))) 
dataNameList 
dataValueList) 

) 

) 


(defobfun  (saveUmpireData  Umpire)  () 
(with-open-file  (rstream  umpireDataFile 
idirection  : output 
:if-exists  :supersede 
:if-does-not-exisi  :create) 
(prinl  dataValueList  rstream) 
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) 

) 


modified  function 


(defobfun  (setUpForReplication  Umpire)  () 

(let  ((crime  (getCurrentTime)) 

) 

(displayNotHostilitiesFlag) 

(selq  hostilitiesp  nil) 


9«SI1t911l91111191911191199999999999999999999999l999999999999999999999999999999999«99t9 


;;  set  up  globalLambdaSigmaValue  and  schedule  event  in 
exponentially  separate  rime  to  i.pdate  its  value 


9999999, 1  I‘  919919199'1999t999l919919999999999999>99999*«)»M9999999999199999199999991999 

(setq  globalLambdaSigmaValue  (♦  globaiSigma  (normal Draw))) 
(addEvent  (make-coast-event 

:time  (+  cTinie  (exponentialDraw  globalMeanlnterval)) 
.object  (self) 

:proced’jre  'updateGloba’LambdaSigma Value 
:data  nil 

rupdaicList  nil)) 

999199999»91919111919191f  91M99999999199999999999f9l9999999t99919199919199999999999999999 

(begiiiFOMdeteciionOpponuniiies  cTime) 

(bcginFalseAlarms  cTime) 

(addEvent  (make-coast-event 

‘.time  (-»-  prehostilitiesDuration  cTime) 

:object  (selO 

rprocedure  'processBeginHostilities 
;data  nil 
:updatcList  nil)) 

(addEvent  (make-coasi-evcni 

:rime  (+  prehostilitiesDuration  hostilitiesDuration 

cTime) 

:object  (self) 
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:procedure  'processEndOfReplication 
:data  nil 
;updateList  nil)) 


) 


) 


(defobfun  (processBeginHosdlities  Umpire)  (cTime  data) 
(setq  hostilitiesp  t) 

(displayHosiiliiiesFlag)) 


(defobfun  (processEndOfReplication  Umpire)  (cTime  data) 
(setq  hostilitiesp  nil) 

(endReplication) 

(displayNotHostilitiesFlag) 

) 


new  function  added  to  update  global  lambda-sigma 
value  in  exponentially  separate  time  interval 


(defobfun  (updateGlobalLambdaSigmaValue  Umpire)  (cTime  data) 
(setq  globalLambdaSigmaValue  (*  globalSigma  (normalDraw))) 
(addEvent  (make-coast-event 

:time  (-^  cTime  (exponential Draw  globalMeanlnterval)) 
object  (selO 

:procedure  'updateGlobalLambdaSigmaValue 
;data  nil 

:updateList  nil)) 

) 


(defobfun  (beginFOMdetectionOpportunities  Umpire)  (cTime) 

(dolist  (sideStnict  (ask  objectManager  forceStructures)) 

(dolist  (searcher  (getf  (forceStruct-instanceList  sideStruct)  :FOMsearchers)) 
(dolist  (target  (getf  (forceStruct-instanceList  sideStruct)  :FOMtargets)) 
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(setNextDetectionOpportunity  ctime  (list  searcher  target)))))) 


target's  radiated  frequency  changes  depending 
which  environment  region  the  target  is  in; 
if  searcher  is  engaging  the  detection  will  rescheduled  with 
exponential  draw  with  target's  detectable  interval  as  mean 


(defobfun  (setNextDetectionOpportunity  Umpire)(ctime  platforms) 
(let*  ((searcher  (car  platforms)) 

(searcherStaie  (ask  searcher  platformState)) 
garget  (cadr  platforms)) 

(targetState  (ask  target  platformState)) 

(rng  (car  (getRngBrg  (platformStatcParms-lat  searcherState) 
(platformStateParms-lng  searcherState) 
(platformStateParms-lat  targetState) 
(platformStateParms-lng  targetState) 
♦greatCircle*))) 

(maxDetectionRange 
(if  (equalp  'submarineTarget 

(platformStateParms-targetType  targetState)) 
maxSubDetecdonRange 
maxSurfaccDetectionRange)) 
(separation  (-  mg  maxDetectionRange)) 


) 


(targetAcoustics  (ask  target  platformAcoustics)) 
(targetEnvRgn  nil) 

(tageiEnvRgnName  nil) 

(envRgnList  (ask  environmentManagcr  envRegionList)) 
(SF  nil)  ;target's  most  detectable  radiated  frequency 
(mean  Detec  tlnterval 

(platformAcousticParms-detectableInterval 

targetAcoustics)) 


* « 

;;  determine  target's  most  detectable  radiated  frequecy 

t » 


(dolist  (thisEnvRgn  envRgnList) 
(when  (ask  environmentManagcr 
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(objectInRegionp  target  (ask  thisEnvRgn  geoRegion))) 
(setq  targetEnvRgn  thisEnvRgn) 

(return) 

)) 

(setq  targetEnvRgnName  (ask  targetEnvRgn  envRgnName)) 
(setq  SF  (first  (getf  (ask  target 

(platformAcousticParms-env  freq-list 

platfoi  nvAcoustics)) 

targetEnvRgnName))) 


;;  check  that  both  platforms  are  active 


(when 

(and  (platformStateParms-activep  searcherState) 
(platformStateParms-activep  targetState)) 


;;  if  searcher  is  engaged  schedule  next  detection  opportunity 
;;  based  on  earlist  possible  time  these  plafforms  can  come 
;;  within  range  or  on  mean  detection  interval  whichever  is 
;;  longer;  now  engaged  target  still  detectable 
::  by  non-engaged  searacher 


(when  (platformStateParms-engagcdp  searcherState) 

(if  (<  ()  (•  separation  (*  maxClosingSpeed 

meanDetcctIntcrval))) 

(addEvent  (make-coast-event 

:time  (+  cTime  (max  meanDctectInterval 

(/  separation  maxClosingSpeed))) 

:object  (selO 

.procedure  'setNextDctcctionOpportunity 
:data  platforms 
;updatcList  platforms) 

) 

(addEvent  (make-coast-event 

:time  (-^  cTime  (exponentialDraw 
meanDctectInterval)) 

robject  (selQ 

rprocedure  'setNextDctcctionOpportunity 
rdata  platforms 
:updateList  platforms) 
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) 

(retum-from  setNextDetectionOpportunity  nil) 

) 

#1 

depending  upon  range;  begin  random  glimpsing, 
or  calculate  earliest  possible  time  these  platforms  can  come  within  range 
(based  on  maxDectionRange  MaxClosingSpeed  and  meanDetectInterval) 
I# 

(if  (<  0  (-  separation  (♦  maxClosingSpeed 

meanDetectInterval))) 

(addEvent  (make-coast-event 

:time  (■»-  cTime  (max  meanDetectInterval 

(/  separation  maxClosingSpeed))) 

:object  (selQ 

:proccdure  'setNextDetectionOpportunity 
:data  platforms 
:updateList  platforms) 

) 

(addEvent  (make-coast-event 

;time  (+  cTime  (exponential Draw 
meanDetectInterval)) 

:object  (selQ 

;procedure  'processDetectionOpportunity 
:data  platforms 
:updateList  platforms) 

) 

) 

) 

) 

) 


(dcfobfun  (processDetectionOpportunity  Umpire)(ctimc  platforms) 
(let*  ((searcher  (car  platforms)) 

(target  (cadr  platforms)) 

(searcherState  (ask  searcher  platformState)) 

(targetState  (ask  target  platformState)) 
(maxDetectionRange 
(if  (equalp  'submarineTarget 
(plaiformStateParms-targetType  targetState)) 

maxSubDetectionRange  maxSuifaceDetectionRange)) 
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(targetAcoustics  (ask  target  platformAcoustics)) 
(meanDetectlnterval 

(platformAcousticParms-detectablelnterval 

targetAcoustics))  ;target's  detetable  interval 

) 

(when  (and  (platformStateParms-activep  searcherState) 
(platformStateParms-activep  targetState)) 


;;  engaged  searcher  searcher  won't  detect  other  targets 
;;  engaged  target  is  detectable  to  other  searchers 


(when  (platformStateParms-engagedp  searcherState) 
(addEvent  (make-coast-event 

:time  (+  cTime  (exponentialDraw  meanDetectlnterval)) 
:object  (selQ 

iprocedure  'processDetectionOpportunity 
;data  platforms 
.•updateList  platforms) 

) 

(retum-from  processDetectionOpportunity  nil) 

) 

(if  (<  (car  (getRngBrg  (platformStateParms-lat  searcherState) 
(platformStateParms-Ing  searcherState) 
(platformStateParms-lat  targetState) 

(platforms  tateParmS'lng  targetState) 
*greatCircle*)) 
maxDetectionRange) 

(if  (getDetectionOppResult  searcher  searcherState 
target  targetState) 

(processTargetDetection  ctime  searcher  searcherState 
target  targetState)) 

) 

(setNexiDetectionOpportunity  ctime  platforms) 

) 

) 

) 
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;  modified  function 


(defobfun  (getDetectionOppResult  Umpire)(searcher  searcherState  target 
targetState) 

(let  ((searcherAcoustics  (ask  searcher  platformAcoustics)) 

(targetAcoustics  (ask  target  platformAcoustics)) 

(envRgnList  (ask  environmentManager  envRegionList)) 

(targetEnvRgn  nil) 

(targetEnvRgnName  nil) 

(SF  nil) 

) 


;;  determine  target's  most  detectable  radiated  frequecy 


(dolist  (thisEnvRgn  envRgnList) 

(when  (ask  environmentManager 

(objectInRegionp  target  (ask  thisEnvRgn  geoRegion))) 
(setq  targetEnvRgn  thisEnvRgn) 

(return) 

)) 

(setq  targetEnvRgnName  (ask  targetEnvRgn  envRgnName)) 
(setq  SF  (first  (getf  (ask  target 

(platformAcousticParms*env-freq*list 

platformAcoustics)) 

targetEnvRgnName))) 


(multiple-value-bind 
(PL  AN  sigma) 

(ask  environmentManager  (gctPL-AN-Sigma  searcher  target)) 
(let*  ((RN  (if  (<  (plaiformStateParms-spd  targetState) 
(platformAcousticParms-s^Thresh 

targetAcoustics)) 

(platformAcousticParms-RNslow  targetAcoustics) 
(platformAcousticParms-RNfast  targetAcoustics))) 
(SN  (if  (<  (platformStateParms-spd  searcherState) 
(platformAcousticParms-spdThrcsh 

searcherAcoustics)) 
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(car  (getf  (platformAcousticParms-SNlist 

searcherAcoustics)  SF)) 
(cadr  (getf  (platformAcousticParms-SNlist 

searcherAcoustics)  SF)))) 

(DI  (car  (getf  (platformAcousticParms-DI-RD-list 

searcherAcoustics)  SF))) 

(RD  (cadr  (getf  (platformAcousticParms-DI-RD-list 

searcherAcoustics)  SF))) 
(FOM  (if  (and  RN  SN  AN  DI  RD) 

(-  RN  (+  (-  (powerSum  SN  AN)  DI)  RD)))) 

(MSE  (if  (and  FOM  PL)(-  FOM  PL))) 


using  global  and  individual  lambda-sigma-value 
to  adjust  MSE 


(SE”(ifMSE” . . . . . 

(-I-  MSE  (-t-  globalLambdaSigmaValue 
(ask  searcher 

(platformAcousticParmS'lambdaSigmaValue 

platformAcoustics)))))) 

) 

(retum-from  getDetectionOppResult 
(if  SE  (if  (<  0  SE)  t  nil)  nil) 

) 

) 

) 

) 

) 


(defobfun  (processTargetDetection  Umpire)  (ctime  searcher  searcherState 

target  targetSiate) 

(let  ((reportDelay 

(ask  searcher  (getReportDelay 

(platformStateParms-targetType  targetState))))) 
(addEvent  (make-coast-event 

:time  (-♦-  cTime  reportDelay) 
robject  (self) 

.•procedure  'transmitDetectionReport 
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;data  (list  searcher  target) 
lupdaieList  (list  searcher  target) 

) 

) 

(when 

(and  (or  (equalp  'ASW-attack  (ask  searcher  mission)) 
(equ^p  'ASW-trail  (ask  searcher  mission))) 

(equalp  'submarineTarget  (ask  target  targetType)) 

) 

(if  (and  hostilitiesp  (equalp  ‘ASW-attack 

(ask  searcher  mission))) 

(processSubSubEngagement  cdme  reportDcIay  searcher 
searcherState  target  targetState) 
(processTrackAndTrail  ctime  reponDelay  searcher 
searcherState  target  targetState 

(ask  searcher  subTrailTime)) 

) 

) 

(when 

(and  (or  (equalp  'ASuW-attack  (ask  searcher  nnission)) 
(equtilp  'ASuW-trail  (ask  searcher  mission))) 
(equsdp  'surfaceTarget  (ask  target  targetType)) 

) 

(if  (and  hostilitiesp  (equalp  'ASuW-attack 
(ask  searcher  mission))) 

(processSubSurfaceEngagement  cdme  reponDelay  searcher 
searcherState  target  targetState) 
(processTrackAndTrail  cdme  reponDelay  searcher 
searcherState  target  targetState 
(ask  searcher  surfaceTrailTime)) 

) 

) 

) 

) 


(defobfun  (processTrackAndTrail  Umpire)  (cTime  delay  searcher  searcherState 
target  targetState  trailTime) 

(when  (and  (ask  searcher  patrolRegion)(ask  target  patrolRegion)) 

(let*  ((iheinterval  (+  (exponential Draw  trailTime)  delay)) 

) 

(ask  searcher  (beginTrail  cdme  theinterval  target)) 
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(ask  target  (beginEngagement  crime  thelnterval  searcher)) 

(addEvent  (make-coast  -event 

:rime  (-^  crime  thelnterval) 

:objeci  (scIO 
rprocedure  'ProcessEndTrail 
rdata  searcher 
:updateList  (list  searcher) 

) 

) 

(addEvent  (make-coast-event 

:rime  (-h  cTime  thelnterval) 
lobject  (selt) 

.•procedure  'processPlatformDisengage 
rdata  target 
rupdateList  (list  target) 

) 

) 

(recordMOEevent  (make-MOEeventStruct 
reventlD  rtrailingPeriods 
robjectl  (ask  searcher  MOEeventObjectIndex) 
:object2  (ask  target  MOEeventObjectIndex) 
rdata  (dround  thelnterval  3))) 

) 

) 

) 


(defobfun  (processSubSubEngagement  Umpire)  (crime  delay  searcher 

searcherState  target  targetState) 

(let*  ((targetClass  (ask  target  (platformStateParms-classType 

platformStatc))) 

(pkList  (ask  searcher  (geif  SubmarinePkAssocList 

targetClass))) 

(seaicheiPK  (car  pkList)) 

(targetPK  (cadr  pkList)) 

(subSubEngagementTime  (ask  searcher 

subSubEngagementTime)) 

(thelnterval  (+  (exponenrialDraw  subSubEngagementTime) 
delay)) 

) 

(if  (not  (and  (numberp  searcherPKXnumberp  targetPK))) 
(retum-from  processSubSubEngagement  nil)) 
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(if  (ask  searcher  (or  (<  currentASWt  ASWtPE) 

(<  currenlClMs  CtMsPE))) 

(retum-from  processSubSubEngagcment  nil)) 

(setq  searcherPK  (ask  target  (expendCounterMeasures 

searcherPK))) 

(setq  targetPK  (ask  target  (expendASWtorpcdos  targetPK))) 

(setq  targetPK  (ask  searcher  (expendCounterMeasures 

targetPK))) 

(setq  searcherPK  (ask  searcher  (expend AS Wtorpedos 

searcherPK))) 

(let  ((target-killed-p  (if  (>  searcherPK  (random  1 .0))  t  nil)) 
(searcher-killed-p  (if  (>  targetPK  (random  1.0))  t  nil)) 

) 

(ask  searcher  (beginEngagement  ctime  theinterval  target)) 
(ask  target  (beginEngagement  ctime  theinterval  searcher)) 
(if  target-killed-p 
(addEvent  (make-coast-event 

:time  (-h  ctime  theinterval) 

:object  (selO 

:procedure  'processPlatformKill 
;data  target 
:updaieList  (list  target) 

) 

) 

(addEvent  (make-coast-event 

:time  (-♦-  ctime  theinterval) 

:object  (selO 

:procedure  'processPlatformDisengage 
:data  target 
:updateList  (list  target) 

) 

) 

) 

(if  searcher-killed-p 
(addEvent  (make-coast-event 

:ume  {+  ctime  theinterval) 
robject  (selO 

:procedure  'processPlatformKill 
:data  searcher 
:updateList  (list  searcher) 
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) 


) 


y  ' 

) 

(addEvent  (make-coast-event 

:time  (-1-  crime  thelnterval) 

:object  (selO 

:procedure  'processPlatfonmDisengage 
:data  searcher 
:updateList  (list  searcher) 

) 

) 

) 


(defobfun  (processSubSurfaceEngagement  Umpire)  (crime  delay 

searcher  searcherState  target  targetState) 

(let  ((ASuWweaponsRange  (ask  searcher  ASuWweaponsRange))) 

(when 

(>  ASuWweaponsRange  (car  (getRngBrg  (platformStateParms-lat 
searcherState) 

(platformStateParms-lng  searcherState) 
(platformStateParms-lat  targetState) 
(platformStateParms-lng  targetState) 
♦greatCircle*))) 

(ask  searcher  (resumePlan  ctime  interceptLegCount)) 

(if  (>  (ask  searcher  currentASuWt)  0) 

(processSubSurfaceAttack  ctime  searcher  searcherState 

target  targetState)) 

(retum-from  processSubSurfaceEngagement) 

) 

(multiple-value-bind 

(CPAdistance  timeToCPA  InterceptCourse) 

(calculateCPA  searcherState  targetState  (ask  searcher 

interceptSpeed)) 

(if  (and  InterceptCourse 

(<  timeToCPA  (ask  searcher  maxASuWInterceptTime)) 

(<  CPAdistance  ASuWweaponsRange) 

) 

(let  ((cpa  (getlatlng  (platformStateParms-lat  searcherState) 
(platformStateParms-lng  searcherState) 
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(*  timeToCPA  (ask  searcher  interceptSpeed)) 
IniercepiCourse 

(platformStateParms-nav  searcherState)))) 

(unless  gib.batchp 


) 


) 

) 

) 

) 


(drawTrack  (list  (list  (platformStateParms-lat 

searcberState) 

(platformStatePanns-lng  searcherState) 

)  cpa)) 

(sleep  .5) 

(drawTrack  (list  (list  (platformStateParms-lat 

searcherState) 

(platformStateParmS'lng  searcherState) 

)  cpa)) 

(ask  searcher  (intemiptPlan  ctime  Interceptcourse 

timeToCPA  target)) 

) 

(ask  searcher  (resumePlan  ctime  interceptLegCount)) 


(defobfun  (processSubSurfaceAttack  Umpire)  (ctime  searcher  searcherState 
target  targetState) 

(let  ((thelnterval  (exponential Draw  (ask  searcher 

subSurfaceEngagementTime)))) 

(ask  searcher  (processSurfaceAttack  target  targetState 

thelnterval)) 

(ask  searcher  (beginEngagement  ctime  thelnterval  uu^get)) 

(ask  target  (beginEngagement  ctime  thelnterval  searcher)) 

(addEvent  (make-coast-event 

:time  (+  ctime  thelnterval) 

:object  (self) 

iprocedure  'processPlatfonnDisengage 
:data  searcher 
rupdateList  (list  searcher))) 

(addEvent  (make-coast-event 

:time  (+  ctime  thelnterval) 
lobject  (selO 
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:procedure  'processPIatformDisengage 
:data  target 
:updateList  (list  target))) 


) 

) 


(defobfun  (transmitDetecdonReport  Umpire)  (ctime  platforms) 

(let*  ((searcher  (car  platforms)) 

(target  (cadr  platforms)) 

(targetState  (ask  target  platformState)) 

(messageType  (case  (platformStatcParms-targetType  targetState) 
(submarineTarget  ‘Sub-Detect-msg) 

(surfaceTarget  ’Surf-Detect-msg))) 

) 

(when  (not  messageType) 

(retuni-from  transmitDetectionReport)) 

(ask  searcher 

(transmit  ctime 

(make-coast-message 
:send-time  crime 
:type  messageType 

xontent  (list  (make-obu-repon 

ctime  targetState)) 

:size  0 

:transmission-path  (list  searcher) 
;transmis$ion-count  0))) 

) 

) 


(defobfun  (processEndTrail  Umpire)  (crime  platform) 
(ask  platform  (endTrail  ctime)) 

) 


(defobfun  (processPIatformDisengage  Umpire)  (ctime  platform) 
(ask  platform  (endEngagement  ctime)) 

) 
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modified  function 


(defobfun  (processPIatformKill  Umpire)  (come  platform) 

(let*  ((platformState  (ask  platform  platformState)) 
(platform-activep  (platformStateParms*activep 

platformState)) 

) 

;  check  if  it  was  killed  already 
(if  (not  platform-activep) 

(retum-from  processPIatformKill  nil) 

) 

(ask  platform  (die)) 

) 

) 


(defobfun  (beginFalseAlarms  Umpire)  (cTime) 

(dolist  (sideStruct  (ask  objectManagcr  forceStructures)) 

(dolist  (searcher  (getf  (forceStnict-instanceList  sideStruct) 

:FOMsearchers)) 

(setNextFalseAlarm  ctime  searcher  ’surfaceTarget) 
(setNextFalseAlaim  ctime  searcher  ’submarineTarget)))) 


9 

;  modified  function 


(defobfun  (setNextFalseAlarm  Umpire)(ctime  searcher  targetType) 

;  check  if  the  searacher  was  killed 

(if  (ask  searcher  (not  (platformStateParms-activep  platformState))) 
(retum-from  processFalseAlarm  nil) 

) 

(if  (>=  0  (ask  searcher  computedFAR)) 

(retum-from  setNextFalseAlarm)) 

(addEvent 
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(make-coast-event 

:time  (+  crime  (ask  searcher  (exponentialDraw 

(/  24  computedFAR)))) 

robject  (self) 

:procedure  'processFalscAIarm 
;data  (list  searcher  targetType) 
rupdateList  (list  searcher) 

) 

) 

) 


modified  function 


(defobfun  (processFalseAlarm  Umpire Xcrime  data) 

(let  ((searcher  (car  data)) 

(targetType  (cadr  data))) 

;  check  if  the  searcher  was  killed 

(if  (ask  searcher  (not  (platformStateParms-activep  platformState))) 
(retum-from  processFalseAlarm  nil) 

) 

(if  (ask  searcher  (platfoimStateParms-engagedp  platformState)) 
(progn 

(setNextFalseAIarm  ctime  searcher  targetType) 

(retum-from  processFalseAlarm  nil))) 

(let  ((reportDelay  (ask  searcher  (getReportDelay  targetType)))) 

(addEvent 

(make-coast-event 

:time  (-»■  cTime  reportDelay) 
lobject  (selO 

iprocedure  'transmitFalseAlarmReport 
:data  data 
.‘updateList  (list  searcher) 

) 

) 

) 

(setNextFalseAIarm  ;time  searcher  targetType) 
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) 


(defobfun  (transmitFalseAlarmReport  Umpire)  (ctime  data) 

(let*  ((searcher  (car  data)) 

(targetType  (cadr  data)) 

(uncertainty  #l(if  (equalp  'submarineTarget  targetType) 
maxSubDetecdonRange 
maxSurfaceDetecdonRange)  SCLI#  30) 
(platformState  (ask  searcher  PlatformState)) 

(lat  (platformStateParms-lat  platformState)) 

(Ing  (platformStateParmS'lng  platformState)) 
(messageType  (case  targetType 

(submarineTarget  'Sub-Detect-msg) 

(surfaceTarget  ’Surf-Detect-msg))) 

) 

(ask  searcher 
(transmit  ctime 
(make-coast-message 
:send-time  ctime 
:type  messageType 

xontent  (list  (make-false-obu-report 

lat  Ing  uncertainty  cdme  targetType)) 

:size  0 

:transmission-path  (list  searcher) 

:transmission-count  0))) 

) 

) 
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;;  CHANGE  LOG: 

;;  new  data  fields  and  input  dialogs  added  for  individual 
;;  lambda-sigma  process; 

;;  individualSigma 
;;  lambdaSigmaValue 
;;  meaninterval 

;;  glimpse  interval  deleted  now  using  target's  detectable  interval 
;;  to  schedule  next  detection  event: 

;;  detectableinterval 

;;  editEnv-Freq  added  for  inputting  environment/frequency  pairs 
;;  of  a  acousticplatform 

;;  setForStart  modified 

;;  initializes  env-freq-list.  lambda-sigma  process  parameters, 

;;  target's  detectable  interval  and  schedules  event  for  updating 
;;  individual  lambda-sigma  value 

;;  updateLambdaSigma Value  added 

;;  update  lambda-sigma  value  and  schedule  event  for  next  update 

(setq  acousticPlatform  (kindof  generalPlatform)) 


modified  function 


(defobfun  (exist  acousticplatform)  (init-list) 

(usual-exist  init-list) 

(have  'platformAcoustics  (make-platformAcousticParms)) 
(have  'spdThresh  12.5) 

(have  'RNslow  135) 

(have  'RNfast  150) 

(have  'sigFreq  nil) 
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(have  'individualSigma  0) 

(have  'lambdaSigmaValue  0) 

(have  'meanlnterval  1(XX))  ;(=l/lainbda) 

(have  'detectablelnterval  1 )  ;platform's  detectable  interval 


(have  'SNlist  ()) 

(have 'DI-RD-list  ()) 

(push  'SNlist  savedVarNatnes) 
(push  'DI-RD-list  saved VarNames) 


(have  'env-freq-list  ()) 

(push  'env-freq-list  saved  VarNames) 


(have  'dataNameList  (append  dataNameList 
'(:LambdaSigmaProcess 

(IndividualSigma  meanlnterval  detectablelnterval) 
:acousticSig 

( spdThresh  RNslow  RNfast)))) 

(have  'dataValueList  (append  dataValueList 
'  (:LambdaSigmaProcess 

(.individualSigma  .meanlnterval  .detectablelnterval) 
:acousticSig 

( .spdThresh  .RNslow  .RNfast)))) 

(have  'dataTypeList  (append  dataTypeList 

'(:LambdaSigmaProcess  (data  data  data) 

:acousticSig  (data  data  data)))) 

(have  'dataTexiList  (append  dataTextList 
'(:LambdaSigmaProcess 

("Individual  Sigma"  "Mean  Interval(  1/Lambda)  " 

"Detectable  Int. ") 

:acousticSig  ("Speed  TTiresh." 

"Rad  Noise-slow"  "Rad  Noise-fast" 

)))) 

(have  'dataTemplateText  (append  dataTemplateTexi 
'(iLambdaSigmaProcess 
(("  Decibels"  "[0  -  20]") 

("  Hours"  "(0.25  -  unlimited]") 

("  Hours"  "(0.25  -  unlimited]")) 
racousticSig 
( 

("Knots"  "(0-30]") 
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("Decibles"  "[0  -  3001*’) 

("Decibles"  "(0  -  300]"))))) 

(let*  ((rb-designators  '(LambdaSigmaProcess-rb  acousSig-rb 
SN-rb  DI-RD-rb  env-freq*rb  frequency-rb)) 
(rb-values  (list  dummyRBdata  duminyRBdata  dummyRBdata 

dummyRBdata  dummyRBdata  sigFreq)) 
(rb~names  '(dummyRBdata  dummyRBdata  dummyRBdata 

dummyRBdata  dummyRBdata  sigFreq)) 

) 

(update-rb~lists  rb-designators  rb-values  rb-names) 


(defobfun  (editSN  acousticPlatform)  () 

(setq  SNIist 

(inputPairs  (getsearchFreqList  (forceStruct-opponentList 

sideStruct))  SNIist 

(list  "self  noise  by  frequency"  "freq"  "slow  spd"  "fast  spd "))) 


(defobfun  (editDI-RD  acousticPlatform)  () 

(setq  Dl-RD-list 

(inputPairs  (getsearchFreqList  (forceStruct-opponentList 
sideStruct))  DI-RD-list  (list  "DI  and  RD  by  frequency" 


"freq"  "Dl"  "RD"))) 


new  function  for  inputting  environment  frequency  pair 
of  a  acoustic  platform 


(defobfun  (editEnv-Freq  acousticPlatform)  () 

(setq  env-freq-list 

(inputPairs  (ask  environmeniManager  (getEnvRgnNames)) 
env-freq-list  (list  "frequency  in  environments" 

"environment"  "frequency . ')  t)) 
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) 


modified  function 


(defobfun  (setForStart  acousticPlatform)  () 

(usual-selForStart) 

(self  (platformAcousticParms-sigFreq  platformAcoustics) 

sigFreq 

(platfomiAcousticPanns-spdThresh  platformAcoustics) 

spdThresh 

(platformAcousticParms-RNslow  platformAcoustics) 

RNslow 

(platformAcousticParms-RNfast  platformAcousiics) 

RNfast 

(platformAcousticParms-SNlist  platformAcoustics) 

SNlist 

(platformAcousticParms-DI-RD-list  platformAcoustics) 

DI-RD-list 

(platformAcousticParms-env-freq-list  platformAcoustics) 

env-freq-list 

(platformAcousiicParms-individualSigma  platformAcoustics) 

individualSigma 

(platformAcousticParms-lambdaSigmaValue 

platformAcoustics) 

(♦  individualSigma  (normalDraw)) 

(platformAcousticParms-meanlnterval  platformAcoustics) 

meaninterval 

(platformAcousticParms-detectablelnterval 

platformAcoustics) 

detec  tablelnterval 

) 

;  schedule  next  update  event  of  lambda-sigma  value 
(addEvent  (make-coast-event 

:time  (+  (getCurrenttime)  (exponentialDraw  (platformAcousticParms- 
meanlnterval  platformAcoustics))) 

:object  (selQ 

rprocedure  updateLambdaSigmaValue 
;data  nil 
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:updateList  nil) 

) 


new  function 


(defobfun  (updateLambdaSigmaValue  acousticPlatform)  (crime  data) 

(setf  (platformAcousticParms-lambdaSigmaValue 

platformAcoustics) 

(*  (platfonnAcousticParms-individualSigma 
platformAcoustics)  (normalDraw))) 

(addEvent  (make-coast-event 

:time  (+  (getCurrenttime)  (exponential Draw  (platformAcousticParms- 
meanlnterval  platformAcoustics))) 

:object  (self) 

:procedure  'updateLambdaSigmaValue 
:data  nil 
lupdateList  nil) 

) 

) 
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CHANGE  LOG: 

makeMPAassignments  modified 
the  selection  of  next  cue  is  the  cue  with  the  largest  ratio  of 
(MPA's  time  on  station/SPA  size); 

getSearchValues  modified 
searach Values  now  contain 
searchTime  (MPA's  time  on  station) 
missionCount 
MPA 
SPA  size 


(setq  ASWOC  (kindof  command)) 


(defobfun  (exist  ASWOC)  (init-list) 

(usual-exist  init-list) 

(have  'allocationinterval  2) 

(have  'allocationStartTime 
(have  'MPAareaSearchMax 
(have  'minTstalnitial 
(have  'minTstaDiversion  2) 

(have  'maxSPAsize 
(have  'assumedSearcherNav 
(have  'assumedTargetVel 
(have  'assumedTargetTau 
(have  'allocationinterval 
(have  'operatingRegion) 

(have  ’AreaSearchRegion) 

(have  'currentAssignmentList) 

(have  'dataNameList  (append  dataNameList 
'(:ASWOCallocationPanns 
(allocationStartTime  allocationinterval 
MPAareaSearchMax 

maxSPAsize  minTstalnitial  minTstaDiversion)))) 
(have  'dataValueList  (append  dataValueList 
'(:ASWOCallocationParms 


0) 

0) 

4) 

50(XX)) 

♦rhumbline*) 

10) 

4) 

12) 
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(,allocationStartTime  .allocadonlnterval 
.MPAareaSearchMax,  maxSPAsize  .minTstalnitial 
.minTstaDiversion)))) 

(have  'dataTypeList  (append  dataTypeList 
'(:ASWOCallocationParms 
(data  data  data  data  data  data)))) 

(have  'dataTextList  (append  dataTextList 
'(:ASWOCallocationPanns 

("ALLOC  START  TIME"  "ALLOC  INTERVAL" 
"#MPAs  AREA  SEARCH"  "MAX  SPA  SIZE" 

"MIN  TSTA  ASSIGNED"  "MIN  TSTA 

DIVERTED")))) 

(have  'dataTemplateText  (append  dataTemplateText 
’(:ASWOCallocationPanns 
(("Hours"  "[0  -  72)") 

("Hours"  "[0  -  72)") 

("MPAs"  "[0  -  100]") 

("SQR  NM"  "[0  -  100,000)") 

("Hours"  "[0  -  72]") 

("Hours"  "[0  -  72)") 

)))) 


(let*  ((rb-designators  '(Operating-Region-rb  Area-Search-Region- 

rb  ASWOCallocation-rb)) 

(rb-values  (list  operatingRegion  AreaScarchRegion 
dummyRBdata)) 

(rb-names  '(operatingRegion  AreaScarchRegion 
dummyRBdata)) 

) 

(update-rb-lists  rb-designators  rb-values  rb-names) 


(defobfun  (setForStart  ASWOC)  () 

(usual-setForStart) 

(setq  currentAssignmentList  ()) 

(if  (not  operatingRegion)  (reium-from  setForStart)) 

(addEvent  (make-coast-event 

;time  (-♦-  (getCurrentTime)  allocationStartTime 
(-  0.1  (random  0.2))  allocationinterval) 
:o.)ject  (self) 
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:procedure  'processMPAAIlocation 

:data  (+  (getCurrentTime)  allocationStartTime 

allocationlnterval) 

:updateList  nil) 

) 

) 


(defobfun  (processMPAAIlocation  ASWOC)  (ctime  data) 
(addEvent  (make-coast-event 

:time  (+  data  (-  0.1  (random  0.2))  allocationlnterval) 
:object  (selO 

rprocedure  'processMPAAIlocation 
idata  (+  data  allocationlnterval) 

:updateList  nil) 

) 

(resetCurrentAssignmentList  ctime) 

(let  ((search Values  (getSearch Values 

ctime 

(getUnassignedCues) 

(append  (getAvailableMPAs) 
(getDivertableMPAs)))) 

) 

(makeMPAassignments  ctime  searchValues) 

) 

(resetAreaSearch  ctime) 

) 


(defobfun  (getUnassignedCues  ASWOC)  () 

(let  ((a’;‘  Jues  (getTacticalPicture  'Subpicture  operatingRegion)) 
(unasjignedCues  ())) 

(dolist  (cue  allCues) 

(if  (isOpenCuep  cue)  (push  cue  unassignedCues)) 

) 

(retum-from  getUnassignedCues  unassignedCues) 

) 

) 


99»>t»9>»yMM9»99l11»fty»f»9f»*»*9*9**»»iy*»ff»**»**»tt****M9#9l999999999199999 
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(defobfun  (isOpenCuep  ASWOC)  (cue) 

(let  ((open  t)) 

(dolist  (Assignment  currentAssignmentList) 
(when  open 

(let  ((assignedCue  (cadr  Assignment))) 

(when  (<  (car  (getRngBrg  (car  assignedCue) 
(cadr  assignedCue) 

(car  cue) 

(cadr  cue) 

♦greatCircle*)) 

(+  (sqrt  (/  (caddr  cue)  pi); 

(sqn  (/  (caddr  assignedCue)  pi)))) 
(setq  open  nil) 

) 

) 

) 

) 

(return-from  isOpenCuep  open) 

) 

) 


(defobfun  (resetCurrentAssignmentList  ASWOC)  (crime) 
(let  ((newList  ())) 

(dolist  (assignment  currentAssignmentList) 

(if  (<  crime  (car  assignment)) 

('push  assignment  newLisi))) 
setq  currentAssignmentList  newList) 

) 

) 


(defobfun  (getAvailableMPAs  ASWOC)  () 

(let  ((thisAswoc  (self)) 

(availableMPAlist  nil) ) 

(mapc  #'(lambda  (x) 

(mapc  #'(lambda  (y)  (push  y  availableMPAlist)) 

(ask  X  (getAvailableMPAs  thisAswoc)))) 

(getf  (forceStruct-instanceList  sideStruct)  :ownMPAs)) 
(retum-from  getAvailableMPAs  availableMPAlist) 
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(defobfun  (getDivertableMPAs  ASWOC)  () 

(let  ((thisAswoc  (self)) 

(divertableMPAlist  nil) ) 

(mapc  #'(lanibda  (x) 

(mapc  #'(lainbda  (y)  (push  y  divertableMPAlist)) 

(ask  X  (getDivertableMPAs  thisAswoc)))) 

(getf  (forceStruct-instanceList  sideStnict)  :ownMPAs)) 
(return- from  getDivertableMPAs  divertableMPAlist) 

) 

) 


(defobfun  (getMPAsOnAreaSearch  ASWOC)  () 

(let  ((thisAswoc  (selQ) 

(availableMPAlist  nil) ) 

(mapc  #’(lambda  (x) 

(mapc  #'(lambda  (y)  (push  y  availableMPAlist)) 

(ask  X  (getNff  AsOnAreaSearch  thisAswoc)))) 
(getf  (forceStruct-instanceList  sideStruct)  ;ownMPAs)) 
(retum-from  getMPAsOnAreaSearch  availableMPAlist) 

) 

) 


modified  function 


(defobfun  (getSearchValues  ASWOC)  (crime  tracks  AvailableMPAs) 
(lef**  ((trackCount  (list-length  tracks)) 

(J^ACount  (list-length  AvailableMPAs)) 

(searchValues  (make-array  (list  trackCount  MPACount)))) 
(dotimes 
(j  MP/vCount) 

(let*- ((fvlPA  (nth  j  AvailableMPAs)) 

(missionCount  (MPAstruct-missionCount  MPA)) 
(currentLat  (^^Astruct-lat  MPA)) 

(currentLng  (MPAstnict-lng  MPA)) 
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(flightSpeed  (ask  (MPAstruct-squadron  MPA)  flightSpced)) 
(missionTime  (if  (MPAstruct-searchingp  MPA) 

(-  (J^Astnict-missionEndTime  MPA)  crime) 
(ask  (MPAstruct-squadron  MPA)  missionTime))) 
(timeThresh  (if  (MPAstruct-searchingp  MPA) 
minT  staDi  version 
minTstalnitial)) 

) 

(dorimes 
(i  trackCount) 

Oet*  ((track  (nth  i  tracks)) 

(spaLat  (car  track)) 

(spaLng  (cadr  track)) 

(spaArea  (caddr  track)) 

(SPArange  (car  (getRngBrg  currentLat  currentLng 
spaLat  spaLng 
assumedSearcherNav))) 

(transitTime  (/  SPArange  flightSpced)) 

(RetumRange  (car  (getRngBrg  (MPAstruct-baseLat  MPA) 

(MPAstruct-baseLng  MPA) 
spaLat  SpaLng 
assumedSearcherNav))) 

(retumTime  (/  RetumRange  flightSpced)) 

(searchTime  (-  missionrime  (+  transitTime  retumTime))) 

) 

(if  (<  timeThresh  searchTime) 

(let  ((effArea 

(+  spaArea 
(expandAOUpos 
assumedTargetVel 
assumedTargetTau 
(+  transitTime  (/  searchTime  2))))) 

) 

(if  (<  effArea  maxSPAsize) 

f  f 

searach Values  now  contain 
;;  searchTime  (MPA’s  rime  on  station) 

;;  missionCount 
MPA 
SPA  size 
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(self  (aref  searchValues  i  j) 

(list  searchTime 
missionCount 
MPA 

(list  spaLat  spaLng  effArea) 

)) 

)))))) 

) 

(rctum-fiom  getSearch Values  searchValues) 

) 

) 


modified  fMnction 


(defobfun  (makeMPAassignments  ASWOC)  (ctime  searchValues) 
(let*  ((dimensions  (array-dimensions  searchValues)) 
(trackCount  (car  dimensions)) 

(MPACount  (cadr  dimensions)) 

(assignmentList  ()) 

) 

(do  ((doneCueing  nil) 

(nextCue  nil  nil) 

(MPA  nil)(theCue  nil)(iCue  nilXjCue  nil)) 

(doneCueing  ()) 

(setq  doneCueing  t) 

(dotimcs 
(i  trackCount) 

(dotimes 
(j  MPACount) 

(if  (aref  searchValues  i  j) 

(if  nextCue 


the  selection  of  cue  is  decided  on  the  ratio  of 
MPA's  time  on  station/SPA-size 


f»*#**ft»*»»ft*>*t**#f****»t»»**t*t»»»»»*f»***>**»i»«»f 

(when  (or  (>  (/  (car  (aref  searchValues  i  j))  (caddr  (cadddr 
(aref  searchValues  i  j)))) 
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(/  (car  nextCue)  (caddr  (cadddr  nextCue)))) 

(and  (=  (/  (car  (aref  searchValues  i  j)) 

(caddr  (cadddr  (aref  searchValues  i  j)))) 
(/  (car  nextCue)  (caddr  (cadddr 

nextCue)))) 

(<  (cadr  (aref  searchValues  i  j)) 

(cadr  nextCue))) 

) 

(sctq  MPA  (caddr  (aref  searchValues  i  j))) 

(setq  theCue  (cadddr  (aref  searchValues  i  j))) 

(setq  nextCue  (aref  searchValues  i  j)) 

(setq  iCue  iXsetq  jCue  j) 

) 


(progn 

(setq  MPA  (caddr  (aref  searchValues  i  j))) 

(setq  theCue  (cadddr  (aref  searchValues  i  j))) 

(setq  nextCue  (aref  searchValues  i  j)) 

(setq  iCue  i)(setq  jCue  j)) 

) 

) 

) 

) 

(when  nextCue 

(setq  doneCueing  nil) 

(push  (list  MPA  theCue)  assignmentList) 

(dotimes  (i  trackCount)  (self  (aref  searchValues  i  jCue)  nil)) 
(dotimes  (j  MPACount)  (setf  (aref  searchValues  iCue  j)  nil))) 

) 

(dolist  (assignment  assignmentList) 

(let  ((P^A  (car  assignment)) 

(cue  (cadr  assignment))) 

(setq  currentAssigmentList 

(remove-if  #'(lambda(xXequalp  x  (MPAstruct-keepOutOrder 

MPA))) 


currentAssignmentList)) 

(ask  (MPAstnict'Squadron  MPA)  (dispatchMPAtoCue  ctime  MPA 

cue)) 


(push  (list  (MPAstruct-missionEndTime  MPA)  cue) 

currentAssignmentList) 


)))) 
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(defobfun  (resetAreaScarch  ASWOC)  (ctime) 

(let*  ((AvailableMPAs  (getAvailableMPAs)) 

(MPAsOnAreaSearch  (getMPAsOnArea^arch)) 
(searchRegion  areaSearchRegion) 

) 

(son  availableMPAs  #'(Iambda(x  y) 

(<  (MPAstnict-missionCount  x) 
(MPAstnict-missionCount  y)))) 

(dotimes 

(i  (min  (-  MPAareaSearchMax  (list-length  MPAsOnAreaSearch)) 
(list-length  AvailableMPAs))) 

(let*  ((MPA  (nth  i  availableMPA:,. 

(search Area  (ask  (MPAstruct  ^^uadron  MPA) 

uncueoSearchArea)) 


) 

))) 


(cuePoint  nil) 

(cue  nil) 

(openCuep  nil) 

) 

(dotimes  (i  20) 

(when  (not  openCuep) 

(setq  cuePoint  (getRandomPointInRegion  searchRegion)) 
(setq  cue  (list  (car  cuePoint)  (cadr  cuePoint) 
searchArea)) 

(setq  openCuep  (isOpenCuep  cue)) 

) 

(if  (and  (=  i  19)(not  openCuep))(print  "overlap")) 

) 

(ask  (MPAstruct-squadron  MPA)  (dispatchMPAtoAreaSearch 

ctime  MPA  cue  searchRegion)) 
(push  (list  (MPAstruct-missionEndTime  MPA)  cue) 

currentAssignmentList) 
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CHANGE  LOG: 


new  data  field  added  in  MPAstruct 
engaged 

dispatchMPAtoAreaSearch  and  dispatchMPAtoCue  modified 
schedules  detection  opportunity  for  the  MPA  to  every  possible 
submarine  target  bas^  on  the  inversion  of  detection  rate  to  the  ;; 

setNextFalseAlarm  added  to  schedule  FA  event 

endMPAflight  modified 

reset  MPAstnict-engagedp  flag  to  NIL  after  finish  its  flight 

beginSearch  modified 

schedules  detection  opportunity  for  the  MPAQ  to  every  possible 
submarine  target  based  on  the  inversion  of  detection  rate  to  the  ;; 

setNextDetectionOpportunity  added 
glimpse  rate  detetection  model  used;  now  target's  radiated 
frequency  is  used  to  get  DI,  RD,  NP.  sweepWidth; 
then  mean  glimpse  interval  is  computed 
by  inversing  the  detection  rate(NVW/A); 
excludes  submarine  target  not  inside  search  area  before  MPA 
search  ends  depending  upon  range:  begin  random  glimpsing, 
or  calculate  earliest  possible  time  these  platforms  can  come 
within  range  (based  on  maxDectionRange  MaxClosingSpeed  and 
Glimpseinterval);  engaged  MPA  won't  have  further  detection  on 
other  targets 

proccssDctcctionOpportunity  added 
if  MPA  is  already  engaged  then  it  schedule  next  detection 
opportunity  bas^  on  target's  acoustic  parameters  (i.e.  radiated 
frequency,  DI,  RD,  NP)  to  get  sweepWidth  and  calculate  detection 
rate  of  target,  then  inven  the  detection  rate  to  get  glimpse 
interval;  otherwise  it  will  check  if  the  target  is  detected  by  it 

targetDetectedp  and  targetcoveredp  added 
check  if  the  target  is  inside  the  MPA  search  circle 

setNextFalseAlarm  added 

schedule  false  alarm  event  during  MPA  on  station  time 


target 


target 
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by  making  a  exponential  draw  from  mean  time  between  false 
alarm  (/  24  computedFAR) 

processFalseAlaim  added 

if  the  MPA  is  still  on  station  and  not  being  killed  it  will  transmit 
a  false  target  message,  schedule  next  false  alarm,  expend  buoys 
for  target  classification 

endEngagement  added 
reset  ^ffAstruct•engagedp  flag  to  NIL 

transmitFalseAlarmReport  added 
transmit  the  false  alarm  report  to  its  command 

localizeTarget  modified 

check  if  the  MPA  is  on  station  and  the  target  is  active 
before  doing  anything 

processDetectionEvent  modified 
check  if  MPA  is  on  station  and  the  target  is  active 
before  doing  anything 

prosec  uteTarget  added 

if  time  permit  and  still  carry  enough  weapons  MPA  will 
stay  for  another  run  after  it  finish  prosecution 

loseMPA  modified 

check  if  the  MPA  has  not  being  killed  yet 
getSubsInSearchArea  MPAsquadron  deleted 


(setq  MPAsquadron  (kindof  communicator  detectionReporter)) 


(defobfun  (exist  MPAsquadron)  (init-list) 

(usual-exist  init-list) 

(have  ’commander) 

(have  ’commandObjeci) 

(have  'commandindex  (symbol-value  (getf  init-list 

xommandindex))) 
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(have  ’MPAlist) 

(have  ’subList) 

(have  'SubmarinePkAssocList  ()) 

(push  'SubmarinePkAssocList  savedVarNames) 

(have  'squadronSize  9) 

(have  'nussionTime  1 1 ) 

(have  'tumaroundTime  6) 

(have  'failureTime  100) 

(have  'repairTime  48) 

(have  'prosecutionTime  2) 

(have  'channelsMonitored  16) 

(have  'flightSpeed  200) 

(have  'uncuedSearchArea  1()000) 

(have  'assumedSearcherNav  ♦rhurrioline'*') 

(have  'buoyLoadOut  50) 

(have  'buoysPerSearch  19) 

(have  'buoysPerEngagement  4) 

(have  'torpedoLoadOut  4) 

(have  'torpedosPerEngagement  2) 

(have  'DI-RD-list  ()) 

(push  'DI-RD-list  savedVarNames) 

(have  'dataNameList  (append  dataNameList 
'^IviPAsquadron 

vSMi.adroifw‘»ize  missionTime  tumaroundTime 
;reTitnc  repairTime  prosecutionTime 
chdnr.ft?:,Moiiitoied  flightSpeed  uncuedSearchArea ) 
:M?AStores 

(buoyLoadOut  buoysPerSearch  buoysPerEngagement 
torpedoLoadOut  torpedosPerEngagement)))) 

(have  'dataValueList  (append  dataValueList 
'(:MPAsquadron 

(,squadronSize  .missionTime  .tumaroundTime 
TailureTime  .repairTime  .prosecutionTime 
.channelsMonitored  .flightSpeed 
.uncuedSeaichArea) 

:MPAStores 

(.buoyLoadOut  .buoysPerSearch 
.buoysPerEngagement,  torpedoLoadOut 
.torpedosPerEngagement)))) 

(have  'dataTypeList  (append  dataTypeList 
'(:MPAsquadron 

(data  data  data  data  data  data  data  data  data) 
:MPAStores 
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(data  data  data  data  data)))) 

(have  'dataTextList  (append  dataTextList 
’(:MPAsquadron 

("Squadron  Size"  "Mission  Time"  "Turnaround  Time" 
"Exp  Time  to  Fail"  "Exp  Time  to  Repair"  "Exp 
Prosecution  Time"  "Total  Channels"  "Transit  Speed" 
"Uncued  Search  Area") 

:MPAStores 

("Buoy  Loadout"  "Buoys  /  Search"  "Buoys  / 
Localization"  "Torpedo  Loadout" 

"Torpedo  /  Engagement")))) 

(have  'dataTemplateText  (append  dataTemplateText 
'(:MPAsquadron 

(("Integer"  "(0  -  50]")("Hours"  "[0  -  500]") 

("Hours"  "[0  -  500)")("Hours"  "[0-500]") 

("Hours"  "[0  -  500]"K"Hours"  "[0  -  24]") 

("Integer"  "[0  -  50]")("NM  per  Hr"  "[0  -  24]") 
("NM  squared"  "[0  -  50000]")) 

:MPAStores 

(("Integer"  "(0  -  100]") 

("Integer"  "[0  -  100]") 

("Integer"  "(0  -  100]") 

("Integer"  "(0  -  100]") 

("Integer"  "[0-  100]"))))) 

(self  (getf  dataTextList  :reports) 

("REPORT  TYPE"  "POSITION  ERROR"  "SPEED  ERROR" 
"COURSE  ERROR"  "Exp  Localization  Delay" 

"False  Alarm  Rate"  "Track  Initiation")) 

(setf  (getf  dataTemplateText  :repons) 

'(("POS  or  POS-t-VEL"  "[POS  ,  POS+VEL]") 

("Nautical  Miles"  "(0  -  1000]") 

("Knots"  "[0  -  30]") 

("Degrees"  "(0-180]") 

("Hours"  "(0-100]") 

("FA  per  Hour"  "(0  -  1  ??]") 

("Single  Contact  to  Track"  "0  -  no,  1  -  yes"))) 

(let*  ((rb-dcsignators  '(MPAsquadron-rb  MPAStores-rb 
commander-rb  DI-RD-rb  MPASubPk-rb  )) 

(rb-values  (list  dummyRBdata  dummyRBdata  commander 
dummyRBdata  dummyRBdata )) 

(rb-names  '(dummyRBdata  dummyRBdata  commander 
dummyRBdata  dummyRBdata )) 
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) 

(update-rb- lists  rb-designators  rb- values  rb-names) 

) 

(remove-rbs  '(iocation-rb)) 

(if  (getf  init-list  :isClass)(remove-rbs  '(commander-rb 

fusionCenter-rb))) 

(when  (not  (getf  init-list  :isClass)) 

(remf  C3-connectivities  ’Surf-Detect-msg) 

(remf  C3-connectivities  'Sub-SI-msg) 

(remf  C3-connectiv!ties  ’Surf-SI-msg) 

(update-rb-lists  '(Connect-rb)  (list  C3-connectivities) 

'(C3-connectivities))) 


) 


(defobfun  (editDI-RD  MPAsquadron)  () 

(setq  Dl-RD-list 

(inputPairs  (getsearchFreqList  (forceStruct-opponentList 
sideStruct))  DI-RD-list  (list  "DI  and  RD  by  frequnecy" 

"freq"  "DI  " 

) 


"RD"))) 


(defobfun  (editMPASubPKs  MPAsquadron)  () 

(let  ((classList 

(getClassesFromIndices 
(getf  (forceStruct-objectKeyUst  sideStruct) 

:submarineTargets)))) 

(setq  SuumarinePkAssocList 

(inputPairs  classList  SubmarinePkAssocList 

(list  "pk  vs.  sub  class"  "class”  "pk"  "target  pk"))) 

) 

) 


(defobfun  (setForStart  MPAsquadron)  () 
(usual -setForS  tart) 

(setq  commandObject  (getObjectFromName 
commander 
commandindex 


Ill 


(platformStateParmS'Side  platformState))) 
(have  'subTargetList  (getf  (forceStnict-instanceList  sideStruct) 

isubmarineTargeis)) 

(if  (not  commandObject)(retum-from  setForStart  nil)) 

(setq  MPAlist  nil) 

(dotimes  (i  squad^onSi^o) 

(push  (make-MPAstruct 

:ASWOC  commandObject 

:  squadron  (self) 

rbaseLat  (ask  commandObject 

(platformStateParms-lat  platformstate)) 
ibaseLng  (ask  commandObject 

(platformStateParms-Ing  platformstate)) 

:lai  (ask  commandObject 

(platformStateParms-lat  platformstate)) 

:lng  (ask  commandObject 

(platformStateParms-lng  platformstate)) 
tbuoyCount  buoyLoadOut 
:torpcdoCount  torpedoLoadOut 
:failureTime  (exponentialEhaw  failureTime) 

) 

MPAlist 

) 

) 

) 


modified  function 


(defobfun  (dispatchMPAtoAreaSearch  MPAsquadron)  (ctime  MPA  cue 
arcaSearchRegion) 

(let*  ((targetLat  (car  cue)) 

(targetLng  (cadr  cue)) 

(targetRadius  (sqrt  (/  (caddr  cue)  pi))) 

(targetRange  (car  (getRngBrg  (MPAstruct-lat  MPA) 

(MPAstruct-lng  MPA)  targetLat  targetLng 
assumedSearcherNav))) 

(transitTime  (/  targetRange  flightSpeed)) 

(RetumRange  (car  (getRngBrg  (MPAstnict-baseLai  MPA) 

(MPAstruct-baseLng  MPA) 
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targetLat  targetLng 
assumed  SearcherNav))) 
(retumTime  (/  RetumRange  flightSpeed)) 
(searchStartTime  (+  ctirne  transitTime)) 

(missionEndTime  (+  ctime  missionTime)) 
(searchEndTime  (-  missionEndTime  retumTime)) 
(flightHours  (+  (MPAstnict-totalFlightHours  VffA) 
missionTime)) 

(misssionCount  (+  (MPAstruct-missionCount  MPA)  1)) 

) 

(unless  glb.batchp 

(if  (MPAstruct-drawnMission  MPA) 

(plotCueingPicture 

(nth  0  (MPAstruct-drawnMission  MPA)) 

(nth  1  (MPAstruct-drawnMission  MPA)) 

(nth  2  (MPAstruct-drawnMission  MPA)) 

(nth  3  (MPAstruct-drawnMission  MPA)) 

(nth  4  (MPAstruct-drawnMission  MPA)))) 
(plotCueingPicture 
(MPAstruct-lat  MPA) 

(MPAstruct-lng  MPA) 
targetLat 
targetLng 
targetRadius)) 

(setf  (MPAstruct-drawnMission  MPA)  (list 

(MPAstruct-lat  MPA) 
(MPAstnict-lng  MPA) 
targetLat 
targetLng 
targetRadius) 

) 

(setf  (MPAstruct-searchingp  MPA)  t 

(MPAstruct-searchStartTime  MPA)  searchStartTime 
(MPAstruct-searchEndTime  MPA)  searchEndTime 
(MPAstruct-missionEndTime  MPA)  missionEndTime 
(MPAstruct-lat  MPA)  targetLat 

(MPAstruct-lng  MPA)  targetLng 

(MPAstruct-searchRegion  MPA)  areaSearchRegion 
(MPAstmct-radius  MPA)  targetRadius 

(MPAstnict-totalFlightHours  MPA)  flightHours 
(MPAstruct-missionCount  MPA)  misssionCount 

) 

(setf  (MPAstruct-keepOutOrder  MPA)  (list 
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(MPAstruct-missionEndTime  MPA)  cue)) 
(addEvent  (make-coast-event 

:time  (-(-  crime  transitTime) 

:object  (self) 
rprocedure  'beginSearch 
;data  MPA 

;updateList  subTargetList)) 

(addEvent  (make -coast-event 

:rime  (+  crime  transitTime) 

:object  (self) 

:procedure  'setNextFalseAlarm 
:data  MPA 
:updateList  nil)) 

(addEvent  (make-coast-event 
:time  missionEndTime 
:object  (self) 

:proctdure  'endMPAflight 
:data  MPA 
:updateList  nil)) 

) 

) 


modified  function 


(defobfun  (dispatchMPAtoCue  MPAsquaJron)  (ctime  MPA  target) 
(let*  ((targetLat  (car  target)) 

(targetLng  (cadr  target)) 

(targetRadius  ( sqrt  (/  (caddr  target)  pi))) 


(targetRange 


(transitTime 

(RetumRange 


(retumTime 


(car  (gctRngBrg  (MPAstruct-lat  MPA) 
(MPAstruct-Ing  MPA) 
targetLat  targetLng 
assumedSearcherNav))) 

(/  targetRange  flightSpeed)) 

(car  (getRngBrg  (MPAstruci-baseLat  MPA) 
(MPAstruct-baseLng  MPa) 
targetLat  targetLng 
assumedSearcherNav))) 

(/  RetumRange  flightSpeed)) 
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(searchStartTime  (+  ctime  transitTime)) 

(missionTime  (if  (MPAstnici-searchingp  MPA) 

(-  (MPAstnict-missionEndTime  MPA)  ctime) 
missionTime)) 

(missionEndTime  (+  crime  missionTime)) 

(searchEndTime  (-  missionEndTime  retumTime)) 
(flightHours  (if  (MPAstnici-searchingp  MPA) 
(MPAstruct-totalFlightHours  MPA) 

(-»-  (MPAstruct-totalFlightHours  MPA) 
missionTime))) 

(misssionCount  (if  (MPAstruct-scarchingp  MPA) 
(MPAstruct-missionCount  h^A) 

(+  (MPAstruct-missionCount  MPA)  1))) 

) 

(unless  gib.batchp 

(if  (MPAsiruct-drawnMission  MPA) 

(plotCueitigPicture 

(nth  0  (Mf'Astruct-drawnMission  MPA)) 

(nth  1  (MPAstruct-drawnMission  MPA)) 

(nth  2  (MPAstruct-drawnMission  MPA)) 

(nth  3  (MPAstruct-drawnMission  MPA)) 

(nth  4  (MPAstruct-drawnMission  MPA)))) 
(plotCueingPicture 
(MPAstruct-lat  MPA) 

(MPAstruci-lng  MPA) 
targetLat 
targetLng 
taigetRadius)) 

(self  (MPAstruct-drawnMissim  MPA)  (list 

(MPAstruci-lat  MPA) 
(MPAstnict-lng  MPA) 
targetLat 
targetLng 
targetRadius) 

) 

(seif  (MPAstnici-searchingp  MPA)  t 

(MPAstruct-searchStartTime  MPA)  searchStartTime 
(MPAstruct-searchEndTime  MPA)  searchEndTime 
(MPAstruct-missionEndTime  MPA)  missionEndTime 
(MPAstruct-lat  MPA)  targp.tLat 

(MPAstruct-lng  MPA)  targetLng 

(MPAstruct-radius  MPA)  targetRadius 

(MPAstruct-searchRegion  MPA)  nil 
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(MPAstnici-totalFlightHours  MPA)  flightHours 
(MPAstruct-missionCount  MPA)  misssionCount 

) 

(self  (MPAstruct-keepOutOrder  MPA)  (list  (MPAsmjct-missionEndTime 

MPA)  target)) 

(addEveni  (make-coast-event 

:time  (+  ctime  transitTime) 

:object  (selO 
:procedure  'beginSearch 
:data  MPA 

:updateL’st  subTargetList)) 

(addEvent  (makv'-coast-event 

:time  (+  ctime  transitTime) 

:object  (selO 

:procedure  'setNextFalseAlarm 
:data  MPA 
:updateList  nil)) 

(addEvent  (makc-coast-event 
:time  missionEndTime 
:object  (selO 
:procedure  'endMPAflight 
:data  MPA 
;updateList  nil)) 

) 

) 


(defobfun  (endMPAflight  MPAsquadron)  (ctime  MPA) 

(if  (or  (not  (MPAstruct-searchingp  MPA)) 

(MPAstruct-killedp  MPA)) 

(retum-from  endMPAflight  nil)) 

(unless  glb.batchp 

(if  (MPAstnict-drawnMission  MPA) 

(plotCueingPicture 

(nth  0  (MPAstruct-drawnMission  MPA)) 

(nth  1  (MPAstruct-drawnMission  MPA)) 

(nth  2  (MPAstruct-drawnMission  MPA)) 

(nth  3  (MPAstruct-drawnMission  MPA)) 

(nth  4  (MPAstruct-drawnMission  MPA))))) 

(self  (MPAstruct-lat  MPA)  (MPAstruct-baselat  MPA) 

(MPAstruct-lng  MPA)(MPAstruct-baseLng  MPA) 
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(MPAstruct-readyp 

(MPAstruct-searchingp 

(MPAstruct-searchStartTime 

(MPAstruct-searchEndTime 

(MPAstruct-missionEndTime 

(MPAstruct-radius 

(MPAstruct-searchRegion 

(MPAstruct-drawnMission 

(MPAstruct-keepOutOrder 

(MPAstnict-nextConiaci 

(MPAstnjct-engagedp 


MPA)  nil 
MPA)  nil 
MPA)  nil 

MPA)  nil 
MPA)  nil 
MPA)  nil 
MPA)  nil 
MPA)  nil 
MPA)  nil 
MPA)  nil 

MPA)  nil 


) 

(if  (>  (MPAstruct-totalFlightHours  MPAy 
(MPAstnict-failureTime  MPA)) 

(addEveni  (make-coast-event 

:time  (+  (getCurrentTime)  tumaroundTime 
(exponentialE>raw  repaiiTime)) 
:object  (selO 
;procediire  'repairMPA 
:data  MPA 


:updateList  nil)) 

(addEvent  (make-coast-event 

:time  (+  (getCurrentTime)  tumaroundTime) 
:object  (self) 

:procedure  'readyMPA 
;data  MPA 


:updatcList  nil)) 

) 

) 


(defobfun  (repairMPA  MPAsquadron)  (ctime  MPA) 
(setf  (MPAstnict-failureTime  MPA) 

(4-  (MPAstruct-totalFlightHours  MPA) 
(exponentialDraw  failureTime)) 

) 

(readyMPA  ctime  MPA) 

) 


(defobfun  (readyMPA  MPAsquadron)  (ctime  MPA) 
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(self  (MPAstruct-readyp  MPA)  t 

(MPAstnict-buoyCount  MPA)  buoyLoadOut 
(MPA.stnict-torpedoCount  MPA)  toipedoLoadOut 

) 


) 


(defobfun  (getAvailableMPAs  MPAsquadron)  (thisAswoc) 

(let  ((availableMPAs  ()) 

) 

(dolist  (x  MPAlist) 

(if  (and  (not  (MPAstruct-killedp  x)) 

(equalp  thisAswoc  (MPAstruct-ASWOC  x)) 
(MPAstruct-readyp  x) 

(not  (MPAstruct-searchingp  x))) 

(push  X  availableMPAs))) 

(retum-from  getAvailableMPAs  availableMPAs) 

) 

) 


(defobfun  (getDivertableMPAs  MPAsquadron)  (thisAswoc) 
(let  ((divertab)eMPAs  ()) 

) 

(dolist  (x  MPAlist) 

(if  (and  (not  (MPAstruct-killedp  x)) 

(equalp  thisAswoc  (MPAstruct-ASWOC  x)) 
(MPAstruct-searchRegion  x) 

(>=  (MPAstruct-buoyCount  x)  (+ 
buoysPerEngagement))) 

(push  X  divertableMPAs))) 

(retum-from  getDivertableh^As  divertableMPAs) 

) 

) 


(defobfun  (getMPAsOnAreaSearch  MPAsquadron)  (thisAswoc) 
(let  ((availableMPAs  ()) 

) 

(dolist  (x  MPAlist) 


buoysPerSearch 
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(if  (and  (not  (MPAstnict-killedp  x)) 

(cqualp  thisAswoc  (MPAstruct-ASWOC  x)) 
(MPAstruct-searchRegion  x)) 

(push  x  availableMPAs))) 

(reium-from  getMPAsOnAreaSearch  availableMPAs) 

) 

) 


(defobfun  (plotSelf  MPAsquadron)  () 

(unless  gib.batchp 

(if  (and  (platFormStatePanns-activep  platFormState) 
(clockIsOnp)) 

(mapc  #'(lambda(MPA) 

(if  (MPAstruct-drawnMission  MPA) 
(plotCueingPicture 

(nth  0  (MPAstruct-drawnMission  MPA)) 
(nth  1  (MPAstruct-drawnMission  MPA)) 
(nth  2  (MPAstruct-drawnMission  MPA)) 
(nth  3  (MPAstruct-drawnMission  MPA)) 
(nth  4  (MPAstruct-drawnMission  MPA))))) 
MPAlist 

) 

) 

) 

) 


modified  function 


(defobfun  (beginSearch  MPAsquadron)  (cdme  MPA) 
(expendBuoys  MPA  buoysPerSearch) 

(dolist  (target  subTargetList) 
(setNextDetectionOpportunity  ctime  (list  MPA  target)) 

) 

) 


9f 


119 


new  function 


(defobfun  (setNextDetectionOpportunity  MPAsquadronXctime  data) 
(let*  ((MPA  (car  data)) 

(target  (cadr  data)) 

(targetState  (ask  target  platformState)) 

(targetAcousdcs  (ask  target  platformAcoustics)) 

) 

(if  (or  (not  (MPAstruct-scarchingp  MPA)) 

(MPAstruct-killedp  MPA) 

(not  (platformStateParms-activep  targetState)) 

) 

(retum-from  setNextDetectionOpportunity  nil)) 


glimpse  rate  detetection  model  used;  now  target's  radiated 
frequency  is  used  to  get  DI,  RD,  NP,  swcepWidth 
then  mean  glimpse  interval  is  computed  by  inversing  the 
detection  rate  (NVW/A) 

I# 

(let*  ((mg  (car  (getRngBrg  (MPAstruct*lat  MPA) 

(MPAstruct-lng  MPA) 

(platformStateParms-lat  targetState) 

(platforms tateParmS'lng  targetState) 
assumedSearcherNav))) 

(separation  (-  mg  (MPAstruct-radius  MPA))) 
(maxClosingSpecd  (ask  umpire  maxClosingSpeed)); 

(RN  (if  (<  (platformStateParms-spd  targetState) 

(platformAcousticParms-spdl'hresh  targetAcoustics)) 
(platformAcousticParms-RNslow  targetAcoustics) 
(platformAcousticParms-RNfast  targetAcoustics))) 
(targetEnvRgn  nil) 

(tagetEnvRgnName  nil) 

(envRgnList  (ask  environmentManager  envRegionList)) 

(SF  nil) 

(DI  nil) 

(RD  nil) 

(NP  nil) 
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(swecpWidth  0) 
glimpselnterval 
meanGlimpseInterval 
(searchTime  2) 

) 

(dolist  (thisEnvRgn  EnvRgnList) 

(when  (objectlnRegionp  target  (ask  thisEnvRgn  geoRegion)) 

(setq  targetEnvRgn  thisEnvRgn) 

(return) 

)) 

(when  targetEnvRgn 

(setq  targetEnvRgnName  (ask  targetEnvRgn  envRgnName)) 

(setq  SF  (first  (getf  (ask  target  (platformAcousticParms-env-freq-list 
platformAcoustics))  targetEnvRgnName))) 

(setq  DI  (first  (getf  DI-RD-list  S^)) 

(setq  RD  (second  (getf  DI-RD-list  SF))) 

(setq  NP  (if  (and  DI  RD)  (-  (-»-  RN  Dh  RD))) 

) 

(setq  sweepWidth  (if  NP  (ask  environmentManager  (getSweepWidth  target 
NP))  0)) 

(if  (or  (eql  0  swecpWidth) 

(eql  0  (*  (ask  target  (cond  (track  meanTransitSpeed) 

(patroIRegion  meanPatrolSpeed) 

(t  interceptSpeed))) 

(min  channelsMonitored  buoysPerSearch) 
swecpWidth)) 

(eql  0  (♦  pi  (MPAstruct-radius  MPA)(MPAstruct-radius  MPA)))) 
(retum-from  setNextDetectionOpportunity  nil) 

) 

(setq  meanGlimpseInterval  (/  1  (/  (**  (ask  target  (cond  (track 
meanTransitSpeed) 

(patroIRegion  meanPatrolSpeed) 

(t  interceptSpeed))) 

(min  channelsMonitored  buoysPerSeaich) 
swecpWidth) 

(*  pi  (MPAstruct-radius  MPA)(MPAsiruct-radius 

MPA))))) 

(setq  Glimpselnterval  (exponentialDraw  meanGlimpseInterval)) 

(setq  searchTime  (-  (MPAstruct-scarchEndTime  MPA)  (MPAstruct- 
scarchStartTime  MPA))) 


121 


#1 

exclude  submarine  target  not  inside  search  area  before 
MPA  search  ends 


I# 


(when  (and  (<  separation  (*  searchTime  maxClosingSpeed)) 

(<  (+  cTime  Glimpselnterval)  (MPAstnict-searchEndTime 

MPA)) 


) 


#1 

check  if  MPA  is  engaged 
I# 

(when  (MPAstruct-engagedp  MPA) 

(if  (<  0  (-  separation  (*  maxClosingSpeed  Glimpselnterval))) 
(addEvent  (make-coast-event 

:time  (+  cTime  (max  Glimpselnterval 

(/  separation  maxClosingSpeed))) 

;object  (selO 

:procedure  'setNextDetectionOpportunity 
:data  data 

:updateList  (list  (selO  target)) 

) 

(addEvent  (make-coast-event 

:time  (+  cTime  Glimpselnterval) 

:object  (selQ 

.‘procedure  'setNextDetectionOpportunity 
:data  data 

:updateList  (list  (selO  target)) 

) 

) 

(retum-from  setNextDetectionOpportunity  nil) 

) 

#1 

depending  upon  range:  begin  random  glimpsing, 
or  calculate  earliest  possible  time  these  platforms  can  come 
within  range  (based  on  maxDectionRange  MaxClosingSpeed  and 
Glimpselnterval) 

I# 

(if  (<  0  (-  separation  (♦  maxClosingSpeed  Glimpselnterval))) 
(addEvent  (make-coast-event 

:time  (+  cTime  (max  Glimpselnterval 

(/  separation  maxClosingSpeed))) 
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robject  (selO 

:procedure  'setNextDetectionOpportunity 
:data  data 

:updateList  (list  (self)  target)) 

) 

(addEvent  (make-coast-event 

:time  (+  cTime  Glimpseinterval) 
lobject  (selO 

:procedure  'processDetecdonOpportunity 
:data  data 

:updateList  (list  (selO  target) 

) 

) 

) 

) 

) 

) 

) 


new  function 


(defobfun  (proccssDeiectionOppoaunity  MPAsquadron)  (ctime  data) 
(let*  ((MPA  (first  data)) 

(target  (second  data)) 

(targetState  (ask  target  platformState)) 
(targetAcoustics  (ask  target  platformAcousdcs))) 

(when  (and  (platformStateParms-acdvep  targetState) 
(MPAstruct-searchingp  MPA) 

) 


(when  (MPAstruct-cngagedp  MPA) 

(let*  ((RN  (if  (<  (platformStateParms-spd  targetState) 

(platformAcousdcParms-spdThresh  targetAcoustics)) 
(pi  atfomiAcousticParms-RN slow  targetAcoustics) 
(platformAcousticParms-RNfast  targetAcousdes))) 
(targetEnvRgn  nil) 

(tagetEnvRgnName  nil) 

(envRgnList  (ask  environmeniManager  envRegionList)) 
(SF  nil) 

(DI  nil) 


123 


(RD  nil) 

(NP  nil) 

(sweepWidth  0) 
glimpselnterval 
meanGlimpseInterval ) 

(dolist  (thisEnvRgn  EnvRgnList) 

(when  (objectInRegionp  target  (ask  thisEnvRgn  geoRegion)) 
(setq  targetEnvRgn  thisEnvRgn) 

(return) 

)) 

(when  targetEnvRgn 

(setq  targetEnvRgnName  (ask  targetEnvRgn  envRgnName)) 
(setq  SF  (first  (getf  (ask  target 

(platfonnAcousticParms-enV'freq-list 

platfoimAcoustics)) 

targetEnvRgnName))) 

(setq  D1  (first  (getf  DI-RD-list  SF))) 

(setq  RD  (second  (getf  DI-RD-list  S^)) 

(setq  NP  (if  (and  D1  RD)  (-  (-1-  RN  Dl)  RD)))) 

(setq  sweepWidth  (if  NP  (ask  environmentManager 

(getSweepWidth  target  NP))  0)) 

(if  (or  (eql  0  sweepWidth) 

(eql  0  (*  (ask  target  (cond  (track  meanTransitSpeed) 
(patrolRegion  nr»eanPatrolSpeed) 

(t  interceptSpeed))) 

(min  channelsMonitored  buoysPerSearch) 
sweepWidth)) 

(eql  0  (*  pi  (MPAstruct-radius  MPA)(MPAstruct-radius 

MPA)))) 

(return -from  processDetectionOpportunity  nil)) 

(setq  meanGlimpseInterval 
(/  1  (/  (*  (ask  target  (cond  (track  meanTransitSpeed) 

(patrolRegion  meanPatrolSpeed) 

(l  interceptSpeed))) 

(min  channelsMonitor^  buoysPerSearch) 
sweepWidth) 

(*  pi  (MPAstruct-radius  MPA)(MPAstruct-radius 

MPA))))) 

(setq  Glimpselnterval  (exponentialDraw 

meanGlimpseInterval)) 


(addEvent  (make-coast-event 

;time  (+  cTime  (exponentialDraw  glimpselnterval)) 
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) 


robject  (selQ 

;procedure  'processDetecrionOpportunity 
:data  data 

;updateList  (if  target  (list  target)) 

)) 


) 

(retum-firom  processDetectionOpportunity  nil) 

) 

(if  (targetDetectedp  MPA  target  targetState  targetAcoustics) 
(processDelecdonevent  ctime  data)) 
(setNextbetectionOpportunity  ctime  data) 

) 


new  function 


(defobfun  (targetDetectedp  MPAsquadron)  (MPA  target  targetState 
targetAcoustics) 

(if  (not  (targetCoveredp  MPA  target  targetState)) 

(retum-froin  targetDetectedp  nil) 

(retum-from  targetDetectedp  t) 

) 

) 


new  function 


(defobfun  (targetCoveredp  MPAsquadronXMPA  target  targetState) 
(retum-from  targetCoveredp 
(<  (car  (getRngBrg  (MPAstruct-lat  MPA) 

(MPAstruct-lng  MPA) 

(platforms  tateParms-lat  targetState) 
(platformStateParms-lng  targetState) 
assumedSearcherNav)) 

(MPAstruct-radius  MPA) 


125 


) 

) 

) 


new  function 


(defobfun  (setNextFalseAlarm  MPAsquadroii)  (cTime  MPA) 
(if  (>=  0  computedFAR) 

(retum-from  setNextFalseAlann)) 

(addEvent 

(make-coast-evenf 

:tiine  (+  crime  (exponentialDraw  (/  24  computedFAR))) 
.object  (self) 

:procedure  'processFalseAlarm 
:data  MPA 
;updateList  (list  (selO) 

) 

) 

) 


new  function 


(defobfun  (processFalseAlarm  MPAsquadron)(ctime  MPA) 

(if  (or  (not  (MPAstruct-searchingp  P4PA)) 

(MPAstruct-killedp  MPA)) 

(retum-ffom  processFalseAlarm  nil) 

) 

(if  (MPAstruct-cngagedp  MPA) 

(progn 

(setNextFalseAlarm  ctime  MPA) 

(retum-from  processFalseAlarm  nil))) 

(let  ((rcponDelay  (getReportDclay  ’submarineTarget))) 

(expcndBuoys  MPA  buoysPerEngagement) ;  expend  buoys  for 

;target  classification 

(addEvent 
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(make-coast-event 

:time  (-^  cTime  reportDelay) 

:object  (self) 

:procedure  'transmitFalseAlannReport 
;data  MPA 
lupdateList  (list  (self)) 

) 

) 

) 

(setNextFalseAlarm  crime  MPA) 

) 


new  function 


(defobfun  (endEngagement  MPAsquadron)  (crime  MPA) 
(setf  (Nff  Astruct-engagedp  MPA)  nil)) 


new  functi  jn 


(defobfun  (transmitFalseAlannReport  MPAsquadron)  (crime  MPA) 
;  check  if  the  MPA  still  on  station  and  not  being  killed 
(if  (or  (not  ('  IPAstruct-searchingp  MPA)) 

(MPAstruci-killedp  MPA)) 

(return-  from  transmitFalseAlannReport  nil) 

) 

(let  ((theReport  (list  (make-false-obu-report 
(MPAstruct-lat  MPA) 

(MPAstruct-lng  MPA) 

(MPAstruct-radius  MPA) 
ctime  'submarineTarget))) 

) 

(transmit  ctime 

(make-coast-message 
:  send- rime  ctime 
.•type  'Sub-Detect-msg 

xontent  theReport 
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) 

) 


:size  0 

•.transmission-path  (list  (self)) 
;transmission-count  0)) 


;;  modified  function 


(defobfun  (ProceseDetectionEvent  MPAsquadron)  (ctime  data) 
(let*  ((MPA  (car  data)) 

(target  (cadr  data)) 

(targetState  (ask  target  platformState)) 

) 

;  check  it  MPA  is  on  station  and  the  target  is  active 
(if  (not  (and  (platformStateParms-activep  targetState) 
(MPAstruct-searchingp  MPA)) 

) 

(retum-from  ProcessDetectionEvent  nil)) 

(setf  (MPAstruct-searchRegion  MPA)  nil) 

(addEvent 
( make-coast-event 

:time  (min  (+  crime  (getReportDelay)) 

(MPAstruct-searchEndTime  MPA)) 

:object  (selO 
;procedure  'localizeTarget 
:data  data 

:updateList  (if  target  (list  target)) 

)) 

) 


;;  modified  function 


(defobfun  (localizeTiirget  MPAsquadron)  (ctime  data) 
(let*  ((Iself  (selO) 
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(MPA  (car  data)) 

(target  (cadr  data)) 

^targetState  (ask  target  platformState)) 
(theP“-ort  (if  target 

(list  (make-obu-report  ctime 
(ask  target  platformState))) 

(list  (make-false-obu-repori 
(MPAstruct-lat  MPA) 
(MPAstruct-lng  MPA) 
(MPAstruct-radius  MPA) 
ctime  *Sub-Detect-msg))))) 


;;  check  if  the  MPA  is  on  station  and  searching 
;;  the  target  is  acdve 


(if  (not  (and  (platformStateParms-activep  targetStatc) 
(MPAsmict-searchingp  MPA)) 

) 

(returi-from  localizeTarget  nil)) 


(transmit  ctime 

(make-coast-message 
:send-time  ctime 
:type  'Sub-Dctect-msg 

rcontent  theRcport 
:size  0 

:transmission-path  (list  (self)) 

:transmission-count  0)) 

(when  (or  (<  (MPAstruct-torpedoCount  MPA)  torpedosPerEngagement) 

(<  (MPAstruct-buoyCount  MPA)  buoysPerErgagement)) 
(endMPAflight  ctime  MPA) 

(retum-from  localizeTarget  nil)) 

(when  (and  (hostilitiesp)  target) 

(expends  uoys  MPA  buoysPcrEngagcment) 

(let  ((engagementTime  (exponentialDraw  prosecutionTiine))) 

(when  (MPAstruct-searchEndTime  MPA) 

(when  (<  (+  cTime  engagementTime)  (MPAstruct-searchEndTime  MPA)) 
(ask  target  (beginEngagement  ctime  engagementTime  IselO) 

(self  (MPAstruct-engagcdp  MPA)  t) 
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(addEvent 

(make-coast-event 

:time  (+  cTime  engagementTime) 
robject  (selQ 
rprocedure  ’prosecuteTarget 
;data  data 

lupdateList  (if  target  (list  target)) 

)) 

) 

) 

) 

) 

) 

) 


modified  function 


(defobfun  (prosecuteTarget  MPAsquadron)  (ctime  data) 

(let*  ((MPA  (car  data)) 

(MPASide  (platformStateParms-side  platformState)) 

(MPAType  (platformStateParms-classType  platformState)) 
(target  (cadr  data)) 

(targetState  (ask  target  platformState)) 

(fargetSide  (platformStateParms-side  targetState)) 
(targetTypc  (platformStateParms-classType  targetState)) 
(pklist  (getf  SubmarinePkAssocList  targetType)) 
(searcherPK  (ask  target  (expendCounterMeasures  (car  pklisi, 
(targetPK  (ask  target  (expendSubSams  (cadr  pklist)))) 
(objeciNumbei  MOEcventObjeciIndex) 

) 

(if  (not  (and  (platformStateParms-activep  targetState) 
(MPAstruct-seaichingp  MPA)) 

) 

(retum-from  prosecuteTarget  nil)) 

(w^en  (not  (and  (numberp  searcherPKXnumberp  targetPK))) 
(ask  target  (endEngagement  ctime)) 

(retum-from  prosecuteTarget  nil)) 

(cxpendTorpedos  MPA  torpedosPerEngagement) 
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(when  (>  targetPK  (random  1.0)) 

(loseNff  A  ctime  MPA  (ask  target  MOEeventObjectIndex))) 


;;  when  MPA  kill  the  target  it  will  start  another  run  if  time 
;;  permit  and  still  enough  weapons  and  buoys 


(when  (>  searcherPK  (random  1.0)) 

(if  (and  (<  crime  (MPAstruct-searchEndTime  MPA)) 

(<  (MPAstruct-torpedoCount  MPA)  torpedosPerEngagement) 
(<  (MPAstruct-buovCount  MPA)  buoysPerSearch) 

) 

(endEngagement  ctime  MPA) 

(endMPAflight  ctime  MPA) 

) 

(ask  arget  (die  objectNumber)) 

(retum-from  prosecuteTarget  nil)) 


(if  (MPAstruct-killedp  MPA)(retum-from  prosecuteTarget  nil)) 

(when  (or  (<  (MPAstruct-torpcdoCount  MPA) 
torpedosPerEngagement) 

(<  (MPAstruct-buoyCount  MPA)  buoysPerSearch)) 
(endMPAflight  crime  MPA) 

(ask  target  (endEngagement  crime)) 

(retum-from  prosecuteTarget  nil)) 

(expendBuoys  MPA  buoysPerEngagement) 

(expendTor^os  MPA  torp^osPerEngagcment) 

(when  (and 

(>  (ask  target  subEscapeMPAprob)  (random  1 ,0)) 

(>  searcherPK  (random  1.0))) 

(ask  target  (die  objectNumber)) 


f  t 

;;  if  rime  permits  MPA  will  stay  for  another  run 

f  9 

«9f99«*ff«99»9999ft999999fff9*l9f999««999f9999 99 999t9« 9*999999999999 

(if  (<  ctime  (MPAstruct-searchEndTime  MPA)) 
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(endEngagement  cdme  MPA) 
(endMPAflight  ctime  MPA) 

) 

(retum-from  prosecuteTarget  nil) 

) 

(ask  target  (endEngagement  crime)) 

) 


(defobfun  (expendBuoys  MPAsquadron)  (MPA  theNumber) 
(setf  (MPAstnict-buoyCount  MPA) 

(max  0  (-  (MPAstnict-buoyCount  MPA)  theNumber))) 
(recordMOEevent  (make-MOl^ventStnict 
:eventID  ibuoyexpense 
robjectl  MO^ventObjectlndex 
:object2  nil 
:data  theNumber)) 

) 


(defobfun  (expendTorpcdos  MPAsquadron)  (MPA  theNumber) 
(setf  (MPAstruct-torpedoCount  NffA) 

(max  0  (-  (MPAstruct-torpedoCount  MPA)  theNumber))) 
(recordMOEevent  (make-MOEeventStnict 
;eveiitID  :MPAtorpedoExpense 
:objectl  MOEeventObjectIndex 
:object2  nil 
•.data  theNumber)) 

) 


modified  function 


(defobfun  (loseMPA  IvlPAsquadron)  (crime  MPA  killer) 
(if  (MPAstnict-killedp  MPA) 

(retum-ffom  loseMPA  nil)) 

(setf  (MPAstnict-killedp  MPA)  t) 
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(endMPAflight  ctime  MPA) 

(let  ((Iside  side)) 

(ask  MOEmonitor  (updateRunningMOE  :MPAkill  Iside))) 
(recordMOEevent  (m^e-MOEeventStruct 
reventID  :attrition 
:objectl  MOEevcntObjectIndex 
:object2  killer))) 
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